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Unter der Voraussetzung positiver Koeffizienten a, der charakteristischen Gleichung läßt sich mit 
Hilfe der Vorzeichenregel von Routh und des Satzes von Orlando zeigen, daß zur Entscheidung der 
Frage nach dem Vorhandensein von konjugiert komplexen Wurzeln mit positivem Realteil von den n 
H ur witz-Determinanten nur jede zweite bekannt sein muß. Insbesondere gilt, daß die Zahl der Vor- 
zeichenwechsel innerhalb der „Haupi“-Folge Hn_,, Hn-3..., worin H; die k-stellige H ur wit z-Deter- 
minante bedeutet und im Falle ungeraden n noch a, anzufügen ist, gleich ist der Zahl der konjugiert kom- 
plexen Wurzelpaare mit positivem Realteil. 


The author discusses the question, if an algebraic equation with positive coefficients a, complex has roots 
with positive real parts. By aid of Routh’s rule of signs and Orlando’s Iheorem, it is possible to show that, 
for the solution of the problem, only euch second ofthe H ur witz determinants must be known. Especially, 


ihe number of variations in sign of the main set H„-ı, Hn-s, - - . is equal to ihe number of pairs of complex 
roots with positive real parts. Hr is the k-th order H ur witz determinant, and, if n is odd, a, must be added 
to the set. 


L’auteur discute la question, si une Equation algebrique avec des coefficients a, positives a des racines 
complewes avec des parties r&elles positives. A l’aide de la rögle de signes de Routh et du theoreöme d’Orlando 
on peut montrer que pour la solution de ce probleme il ne faut connaitre que chaque seconde des n determinantes 
d@Hurwitz. Spezialement le nombre des variations de signes au milieu de la suite ‚„‚principale“ Hy„-ı, 
Hn-3: . .. (0% Hr signifie la determinante #’H urwitz du k-i&me ordre et, sin est impair, a, doit ötre 
ajout£) est &quivalent au nombre des paires de racines complexes avec des parties reelles positives. 


B nacroameli cTaTbe oÖcy:kaeTca BOHPOC O HAJAIMYHH COUPAIKCHHO-KOMINIEKCHEIX KOPHeH 
C HONOKHTENBHOH NeÄCTBUTENAIBHOH YACTbW. B IPeANONOMKeHHH HOJNOMKHTENIBHEIX KO9hPHLMEHTOB 
4, KAPAKTEPHCTUYeCKOTO YPABHEHHA MOIKHO IIOKASATB IIPH IOMOINH mpaBuıa 3HaKkoB Payca u 
TeopeMst OPAaHA0, YTO, AA pemeHua 9TOTO BOIMPOCA, IHIIB KAKIBIH BTOPOH M3 n ONpenesm- 
rerıeü ['ypsuna AomkeH ÖbITb UsBecteH. B MACTHOCTH CIPABAEAIHBO IIONOKEHHEe, YTO YHCAO 


tHepeMeH 3HaKa B IpeNe1ax ”’TIIaBHOH” TOCHENOBATEJIBHOCTH Hy-ı, Hn-3, ..., TAe Hr — k-sHa- 


YHBIH OImpenernureis I'ypBuba, U, B cıyyae HEeYeTHOTO n, IIPHCOBOKYILIAETCH EINE a,, PABHO 
YHCAHy COIPA:KEHHO-KOMINIERKCHBEIX IIAP KOPHeH C IIONOKHTENIBHOH MeÄCTBUTENIBHOH YACTbI. 
1. Voraussetzung und Fragestellung 
Bekanntlich lassen sich die Differentialgleichungen eines linearen Regelungssystems 
stets befriedigen durch den Exponentialansatz: 
I ROHR ABA LER UT CH), 
wobei sich dann die n Eigenwerte A ergeben aus der charakteristischen Gleichung »-ten Grades: 
aA" + a, Am! Ara al Ola iin aneinsy. (2). 
Es bedeutet keine Einschränkung, wenn wir vereinbaren, daß diese Gleichung jedenfalls so 
geschrieben werden soll, daß 
ln ln er u re 6), 


ist. 

Die Eigenwerte A können einmal] reell sein, es kann sich also um eine monolone Eigen- 
bewegung handeln. In diesem Falle wollen wir schreiben 
Oder es treten paarweise konjugiert komplexe Eigenwerte zusammen auf, die auf oszillatorische 
Bewegungen führen, die wir dann ausdrücken können durch Abklingkonstante ö und Kreis- 
frequenz w: E 

ee) 1 = 1 one BR PR Er BE Er Bu Eu u zu AGJR 

Das System ist als stabil anzusehen, wenn sämtliche o- und ö-Werte positiv sind, oder zusam- 
mengefaßt, wenn der Realteil sämtlicher Eigenwerte negativ ist: 


nr nenne tan» (Ö) 
17 
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Das sind im allgemeinen n-Bedingungen, weshalb die daraus abgeleiteten einfacher zu berech- 
nenden Stabilitätskriterien auch immer als n unabhängige Bedingungen in Erscheinung treten. 

Wenn allerdings oszillatorische Eigenbewegungen auftreten, so reduziert sich die Zahl 
der Fragestellungen, weil ja dabei immer zwei Eigenwerte durch dieselbe Abklingkonstante ö 
hinsichtlich der Stabilität gekennzeichnet sind. 

Nun können wir freilich der charakteristischen Gleichung im Allgemeinen nicht ansehen, 
wieviel reelle und wieviel konjugierte komplexe Wurzeln sie hat. Dagegen können wir sofort 
sehen, ob sie eine positiv reelle Wurzel] haben kann. Wenn nämlich sämtliche Koeffizienten der 
charakteristischen Gleichung positiv sind, kann es keine positiven reellen Wurzeln geben (also 
keine negativen o-Werte), da dabei auf der linken Seite lauter positive Summanden auftreten 
würden. Wir können also zuordnen: 


alle a4> Oma > 02 2. Sr 


(Dieser Zusammenhang stellt im übrigen den einfachsten Spezialfall der bereits von Descartes 
aufgestellten Vorzeichenregel dar, nach welcher höchstens soviele positiv reelle Wurzeln in 
einer algebraischen Gleichung auftreten, als die Koeffizientenreihe Vorzeichenwechsel aufweist). 
Da die Forderung positiver a, zugleich eine notwendige Stabilitätsbedingung für mono- 
tone und oszillatorische Eigenbewegungen ist, wie man an der Aufspaltung von (2) in die 
Faktoren: 
VA + aA... ++ = 
a0 (A+ 01) (A+ 0) - At +8) (AH). ee») 


erkennt, die außer A nur positive Größen enthalten, ist dies zugleich der praktisch haupt- 
sächlich interessierende Fall. 
In diesem Falle können aber beiungeradem.n, wo (2) jedenfalls eine reelle Wurzel 
hat, höchstens nur noch 
(m -—1)/2 wen: von + eig MAR ne, 2 (9a) 


instabile oszillatorische Eigenbewegungen auftreten. Bei geradem n lautet die entsprechende 
Zahl nicht einfach n/2, sondern sie verringert sich noch um die Möglichkeit, daß sämtliche 
Wurzeln zu oszillatorisch instabilen Vorgängen gehören. 

Wie nämlich der Vergleich der Glieder vom Grade A*—1in (8) zeigt, ist allgemein 


4u=%(&0+2 26) drassiran dere, Battagzers su. (10). 

a, >0 schließt also bei geradem n das alleinige Auftreten von negaliven ö-Werten aus, wir 
haben also nur 

(n/2 — 1)L "Ballen ZU 2 5 2er [3b 


Fragestellungen. 

Es ergab sich somit die prinzipielle Frage, ob sich nicht unter dieser fast immer erfüllten 
Voraussetzung die Zahl der verbleibenden Stabilitätskriterien entsprechend erniedrigt. Dies 
trifft z. B. zu, bei dem vom Verfasser an anderer Stelle !) aufgestellten Lagekriterium, welches 
besagt, daß die Schnittpunkte der reellen Achse mit der Ortskurve, die man erhält, wenn 
man in dem Polynom der Gl. (2)A=i.2 setzt’ und 2 von 0 bis oo wachsen läßt, abwechselnd 
rechts und links vom Nullpunkt liegen müssen. 

Auch bei den » von Hur witz?) formulierten Stabilitätskriterien, die bekanntlich 
dadurch gegeben sind, daß sämtliche Determinanten der Folge 


IR 4,0 
170 . 

H=a; H, = BR: H,—=10a,0,0. bis A, od, Eier), 
a 0544435 


in welcher jeweils a, — 0 zu setzen ist, wenn entsprechend dem leicht überblickbaren Bildungs- 
gesetz g> n wird, positiv sind: 


H, 208 Del 2,3 .00Nn 00 az 
auch bei diesen Kriterien lassen sich einfache Beispiele dafür angeben, daß für alle a,>0 
nicht alle 7, untersucht werden müssen. So beinhaltet bein = 4 

44,0 
H,;= |a,a,a,| = a,H,— oda, >0 
0a,a;| 


!)L.Cremer,Z. angew. Math. Mech. 25/27 (1947), Heft 5/6. 
®)Hur witz, A., Mathem. Ann, 46 (1895), S. 273. 
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bereits, daß auch 4, >0Oist. Ebenso weist M. Tolle 3) bereits darauf hin, daß bein — 5 mit 
2720 "und, As 20 


bei Voraussetzung positiver a, bereits H,> 0 erfüllt ist. 

Diese Ersparnis an Kriterien interessiert übrigens nicht nur für die Frage der Stabili- 
tät, wenngleich dies praktisch die wichtigste ist, sondern auch für die allgemeinere Frage, 
wieviel instabile oszillatorische Bewegungen auftreten. 

Hierfür hatte bereits Routh) eine Regel angegeben, die an die Descartes sche 
Zeichenregel erinnert. Sie gilt für beliebige charakteristische Gleichungen, also ohne Beschrän- 
kung auf nur oszillatorisch instabile Systeme und besagt, daß die Zahl der Wurzeln mit posi- 
tivem Realteil gleich ist der Zahl der Zeichenwechsel in der Folge seiner n ‚„‚Probefunktionen‘“, 
vor die noch das positive a, an den Anfang tritt. Routh gab für die Errechnung- dieser 
Funktionen einen einfachen Algorithmus an. Für die folgenden Betrachtungen machen wir 
davon Gebrauch, daß jede Routhsche Probefunktion gleich dem Quotienten aus zwei 
benachbarten Hur wit z -Determinanten ist:>) €) ?) 


H,, 2 Hr \ H,„_. 4 
R, Rz H,_ı — (dy; 1a rg His’ Ben Tr ja b DW mE TER (13), 
und außerdem 
der er H; = 4, ae NE ee (13a). 


Im Spezialfall der Stabilität nimmt also die Routh sche Vorzeichenregel die spezielle Form 
an: 
BD WED il, 2 in Kerr ra ner ustnaa (14), 


führt also im allgemeinen auch auf » Stabilitätskriterien. 


2. Behauptung 
Wir behaupten nun, daß die gestellte Frage zu bejahen ist, indem bei positiven Koeffi- 
zienten der charakteristischen Gleichung bereits mit Stabilität gerechnet werden kann, wenn 
nachgewiesen ist, daß entweder nur die geradstelligen oder nur die ungeradstelligen Hur wit z- 
Determinanten positiv sind. Die jeweils benachbarten ungerad- bzw. geradstelligen sind es 


' dann zwangsläufig. Da stets die eine der beiden Determinantenfolgen von niedrigerer Stellen- 


zahl ist, nämlich allgemein die Folge 
An bar De PELEE REIT DEE (KON), 


die bei geradem n mit dem stets positiven H, endet und bei ungeradem » noch durch das stets 
positive a, zu ergänzen ist, wird man praktisch bei der Verwertung dieser Regel immer 
auf diese einfachere Reihe zurückgreifen; wir wollen sie daher im Folgenden als Haupt- 
folge bezeichnen, und die Folge 
a, für n gerade 
Enaa Hess a, für n ungerade 
als Nebenfolge. 

Wir wollen gleich an dieser Stelle darauf hinweisen, — um einer späteren Enttäuschung 
des nur an der Anwendung interessierten Lesers vorzugreifen — daß der Weg zu der am müh- 
samsten zu berechnenden (n — 1) stelligen-Determinante H,„_ı ganz von selbst an allen H- 
Determinanten von niedrigerer Stellenzahl vorbeiführt, mag man unmittelbar an ihre Aus- 
rechnung herangehen, oder den praktisch vorzuziehenden Weg über den Routhschen Algo- 
rithmus einschlagen. 

Die Rechenersparnis bedeutet also jedenfalls keine Halbierung der Rechenarbeit, sondern 
vielleicht nur eine gewisse Einsparung an Schreibarbeit und anzugebenden Resultaten. 

Dies war bereits zu Beginn der hier mitgeteilten Untersuchung klar. Denn aus einem 
wohl zuerst von Orlando ®)ausgesprochenen Satz, welcher besagt, daß: 


n—1,n 
H,ı=an DI (—-4,—%4) Dee... da... (16) 
H,V 


ist, folgt, daß von allen Hur wit z-Determinanten H,„_, als erste verschwindet, wenn man das 
Gebiet oszillatorisch stabiler Vorgänge verläßt. 


3) Tolle, M., Regelung der Kraftmaschinen. 3. Aufl. Berlin 1921, S. 773. 

4)Routh, E. J., Onthe stability of given state of motion. Adams Price ersay (1877); vgl. auch dessen 
Dynamik der Systeme starrer Körper II, Leipzig (1898). 

5)Bader, W., Arch. f. Elektrot. 34 (1940), S. 293. 

6)Bilharz, H.,Z. angew. Math. Mech. 21 (1941), S. 96. 

?)Cremer, H., Z. angew. Math. Mech. 25/27 (1947), Heft 5/6. 

®$)Orlando, Math. . 71, 8. 233. 
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Die Faktoren des Produktes in (16) lassen sich, sofern A, und A, komplexe, aber nicht 
konjugierte komplexe Wurzeln sind, zu der stets positiven Größe 


(Ay — Ar) (A Ar) Au — Ar) (Au — hr) | 
— [(d. + 8)’+ (ou+ @)?] [(öu+ 6)?+ (wu — @r)?] ) 


zusammenfassen. Ebenso führt die Möglichkeit, daß eine Wurzel reell die andere komplex ist 
zu dem positiven Teilfaktor: 


HSIEMARHTR LO) 


(—A,— A) u h)=(lu tb” +: 0.0.0. (166). 
Unter der Voraussetzung (7) werden auch die Faktoren, in denen beide Wurzeln reell sind, 
positiv: 
(A,—h) E30 + al u a erileg ah Ba lbor 
Einzig der Fall, daß A, und A, ein konjugiertes komplexes Wurzelpaar darstellen 
(I, year ir BRRLEAN 


führt auf Teilfaktoren, die im Rahmen der vorliegenden Fragestellung positiv oder negativ sein 
können. H,„_ı ist also abgesehen von stets positiven Faktoren dem Produkt aller vorkommen- 
den Abklingkonstanten proportional. Hieraus folgt einmal, daß beim Verlassen des stabilen 
Gebietes immer A, das Vorzeichen wechselt. Die Bedingung 


HR Anke a aa 


bestimmt also die Umgrenzung des Gebietes oszillatorisch stabiler Bewegungen. Man kommt 
also sicher nicht darum herum, H,„_ı oder die entsprechende Probefunktion R,_ı auszu- 
rechnen. 

Man übersieht aber auch, daß positives H,_ı allein noch keinen Anhaltspunkt für Stabi- 
lität darstellt, sondern nur besagt, daß die Zahl der oszillatorisch instabilen Eigenbewegungen 
gerade ist. 

Das Ziel der vorliegenden Untersuchungen, welche in Kenntnis dieser Zusammenhänge 
in Angriff genommen wurden, lag daher mehr in dem Nachweis, daß die Zahl der notwendigen 
Kriterien sich mit der Zahl der Fragestellungen in Einklang bringen läßt. 


In dieser Hinsicht können wir zunächst feststellen, daß die Zahl der nun noch zu be- -» 


achtenden H-Determinanten der Hauptreihe sich in Einklang mit (9a) und (9b) bei ungeradem n 
auf (n — 1)/2, bei geradem n auf (n/2) — 1 reduziert hat, da hier H, = a, bereits zu unseren 
Voraussetzungen gehört. 

Ferner ist es mehr vom prinzipiellen als vom Standpunkt der Anwendung aus beachtens- 
wert, daß auch die Frage, wieviel oszillatorisch instabile Eigenbewegungen auftreten, aus dem 
Vorzeichen dieser (n — 1)/2 bzw. (n/2 — 1) H-Determinanten abgelesen werden kann und ins- 
besondere ist es erfreulich, daß sich auch diese Schlußfolgerung in die Form einer Zeichenregel 
fassen läßt. Wir werden nämlich zeigen, daß die Zahl der vorhandenen oszil- 
latorisch instabilen Eigenbewegungen, also der konjugiert komplexen 
Wurzelpaare mit positivem Realteilgleich ist der Zahl der Zeichenwechsel 
inder Hauptfolge (oder auch Nebenfolge) der Hurwitzdeterminanten, 
die bei ungeradem n noch durch das stets positive a, zu ergänzen ist. 


3. Beweis 


Um die Beweisführung möglichst anschaulich zu machen, wollen wir uns einer graphi- 
schen Darstellung bedienen, bei welcher wir zunächst die Folge der H-Determinanten durch 
nebeneinander gezeichnete Kreise darstellen, die je nachdem ob sie positiv oder negativ sind, 
auf einer oberen oder unteren Linie eingezeichnet werden. Ebenso wollen wir die zur Haupt- 
folge gehörigen Kreise groß, die zur Nebenfolge gehörigen klein zeichnen. Diese Kette der 

4 d H-Werte beginnt links mit den stets positiven Grö- 

Y rs Pas Ar Fb Ben a,und a,; da vondi j hd d 
+ © ©&&+-—--0, 0,019 ,0, 0,0% 0 13 6 iesen, je nachdem n gerade 
H | A oder ungerade ist, bald die eine, bald die andere 
zur Hauptreihe gehört, sind sie in Bild 1 mit gro- 
ßen und kleinen Kreisen bezeichnet. Später werden 
er wir uns mit einem großen Kreis begnügen, der ledig- 
Anz Rn lich andeuten soll, daß dieser Anfangswert stets po- 
Bild 1. Kriterien-Schema des stabilen Zustands sitiv ist. Rechts endet diese Kette der H-Werte mit 

- H,„, also mit einem kleinen Kreis. 

Ebenso zeichnen wir darunter die Ro ut hschen Probefunktionen (R-Werte) als Kreuze, 
und zwar auch in eine obere oder untere Zeile, je nachdem sie positiv oder negativ sind. Da 
diese aus dem Quotienten zweier benachbarter H-Werte bestehen, ist es sinnvoll, sie unter die 
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Lücke dieser H-Werte zu zeichnen. Nur am linken Rand, wo die Größen der H- und R-Folge 
identisch gleich a, bzw. a, werden, gehören sie sinngemäß direkt untereinander. 

Bild 1 zeigt zunächst den einfachsten Fall, nämlich den der Stabilität, wobei alle H- und 
R-Werte auf der Plus-Linie liegen. Wir können diesen Fall zum Ausgangspunkt machen, weil 
sich jeder Zustand mit noch so vielen oszillatorisch instabilen Eigenbewegungen durch stetige 
Variation der a-Werte, d.h. aber auch der konstruktiven Daten, von ihm aus stetig erreichen 
läßt und zwar auch ohne daß irgendein a-Wert negativ wird. Dabei kann man ferner immer 
so vorgehen, daß nacheinander erst eine, dann zwei, dann drei oszillatorische Instabilitäten 
auftreten, bis die jeweils interessierende Zahl erreicht ist. Der an sich mögliche Fall, daß 
zwei solche Instabilitäten gleichzeitig auftreten, läßt sich immer durch beliebig kleine Variatio- 
nen umgehen, er stellt also einen Grenzfall dar, der nur dann auftreten kann, wenn sich die 
von uns zu beachtenden Nachbarzustände herstellen lassen. 

Nach der erwähnten Vorzeichenregel von Routh lassen sich die jeweils vorhandenen 
oszillatorischen Instabilitäten aus der Zahl doppelter Zeichenwechsel in der R-Kette ablesen. 
Und da in dieser Kette das erste und letzte Glied stets auf der oberen für positive Werte gelten- 
den Zeile liegen, können wir anschaulich sagen, die Zahl der negativen ö-Werte ist gleich der 
Zahl der Täler in der R-Kette, deren Glieder wir deshalb auch durch Linien verbinden wollen. 

Um nun zu zeigen, daß sicherlich nicht jede beliebige Vorzeichenkombination in den 
H-Werten möglich ist, wollen wir uns einmal als erste Abweichung von dem stabilen Zustand 
den Fall des Bildes 2a ansehen, bei welchem an- 
genommen ist, daß irgendein Wert aus der H- 
Nebenfolge als erster negativ wird. Dies würde, 
wie darunter gestrichelt angedeutet ist, zu einem 
R-Tal, also einer oszillatorisch instabilen Eigen- 
bewegung führen, die wir nicht den Werten der 
Hauptfolge ansehen könnten. Diese Änderung 
ist aber mit dem im zweiten Abschnitt er- 
wähnten Satz von Orlando unverträglich. Hier- 
nach muß, wenn sich die Zahl der oszillatorisch Bild 2. Entstehung eines R-Tales 
instabilen Figenbewegungen ändert, immer By alohb mägliabe „ bimmögliche 
H,_ı und somit auch H, das Vorzeichen wech- 
seln. Das heißt aber, jede neue Instabilität, oder anschaulicher gesagt, jedes neue R-Tal, kann 
nur am rechten Ende entstehen (Siehe Bild 2b). 

Von hier aus kann das Tal nach links zu niedrigeren R-Werten weiter wandern. Der ein- 
fachste Fall ist der, daß die einzelnen H-Werte nacheinander von + nach — rücken, wie das in 
dem Bild 3a und 3b gezeigt ist. Es kann aber auch ein H-Wert übersprungen werden, wie das 
in Bild 4a gezeigt ist. Dann entsteht zunächst kein Wandern, sondern eine Verbreiterung des 


ya Or: zo) %) ©) 6) Q [e) 
gi Omen) O (6) o0 i ai Ö Oo 


a) 


+ O 0.0.0 +0 ae ET 
= (@) 1, ©%0 - ©49:.0,10;; b) 
L b) 

+ wer nenn er = LIE + Kr rennen nee In RZ 
a) 7 
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Bild 3. Einzahniges Weiterwandern eines Bild 4. Talverbreiterung um 2 Zähne und 
einzahnigen Tales Raupenbewegung 


Tales, die aber auf diese Weise nur immer um zwei Zähne erfolgen kann. Geht dann an- 
schließend, wie das in Bild 4b gezeigt ist, der übersprungene H-Wert erst von + nach —, so 
stellt sich das ganze wie eine Raupenbewegung dar, bei welcher das Tal unter abwechselnder 
Verbreiterung und Zusammenziehung immer um zwei Zähne nach links rückt. a 
Jedenfalls können wir uns leicht davon überzeugen, daß alle auf diese Weise möglichen 
einzelnen Täler, die immer nur eine ungerade Zähnezahl haben können, und die wir darum kurz 
als ungerade Täler bezeichnen wollen, in Einklang sind mit der oben aufgestellten Zeichenregel 
der Haupt- oder Nebenfolge der H-Determinanten. Dabei sei aber ausdrücklich darauf hinge- 
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wiesen, daß die in der Hauptfolge und die in der Nebenfolge auftretenden Zeichenwechsel be- 
liebig weit auseinander liegen können. 

Das in dem Bild 3 und 4 gezeigte Weiterwandern des ersten Tales nach links hin ist sogar 
notwendig, wenn ein erneuter Zeichenwechsel in den Werten H,„_-ı und AH, zur Entstehung 
eines zweiten R-Tales führen soll. Es müssen ja allgemein zwei Täler immer wenigstens durch 
einen positiven R-Wert, d.h. aber durch ein H-Wertepaar von gleichem Vorzeichen getrennt 
sein. Bild 5 zeigt in diesem Sinne die Entstehung des zweiten R-Tales. 


a 
Pa SELTEN: (ENT OarG7o 2° On Oo0O0o 2 Oo 
- Oo Ö () sw ORAL E 
Bee na ae er Nast ya x 
S 1 
R RT E RIM 
Bild 5. Entstehung eines zweiten Tales Bild 6. Zahnaustausch zwischen ungeraden Tälern 


Sobald zwei Täler vorhanden sind, ist ein grundsätzlich neuer Vorgang möglich, nämlich 
der Zahnaustausch zwischen zwei sich berührenden Tälern. Ein solcher ist in Bild 6 gezeichnet. 
Die Täler sind zunächst (ausgezogene Linien) ihrer Entstehung nach ungerade. Durch das 
Heraufwandern des Wertes H,„_s verliert das dreizahnige Tal einen Zahn und das einzahnige 
gewinnt einen. Man erkennt, daß alle solchen Austauschvorgänge dazu führen müssen, daß die 
Gesamtzahl der Talzähne erhalten bleibt, daß also aus den beiden ungeraden Tälern zwei gerade 
werden müssen. Wir vermerken weiter, daß durch diesen Vorgang nichts an der Zahl der 
Zeichenwechsel in der Haupt- oder Nebenfolge der H-Determinanten geändert wird. 


Wir wollen nun wieder ein Beispiel einfügen, bei welchem wir zunächst eine Verletzung 
der Zeichenregel konstruieren, um zu erkennen, daß wir dabei in einen Widerspruch geraten. 
Wir fragen uns einfach, was hindert, daß ein solches zweizahniges Tal um einen Zahn nach links 
weiter wandert, was ja (siehe Bild 7a!) zur Folge hätte, daß in der Hauptfolge zwei Zeichen- 
wechsel verschwinden, in der Nebenfolge zwei Zeichenwechsel entstehen würden, indem es 

N hierzu notwendig wäre, daß ein Glied der Haupt- 

HH’ TOT I Oi folge und ein Glied der Nebenfolge gleichzeitig 
2 € i und im umgekehrten Sinne ihr Zeichen wech- 
RE rn Armes seln. Bei einem solchen Doppelwechsel müssen 
A ;. BEN a SEI wir uns fragen, ob der Doppelschritt auch in zwei 
\ Einzelschritten möglich ist, Es läßt sich nun 

im vorliegenden Falle leicht zeigen, daß die Ein- 

zelschritte, gleichgültig, in welcher Reihenfolge 

wir sie vornehmen, nicht möglich sind. Bewegen 


+ O-----.- Oo je) o ie 9.0) K9ro wir zuerst den großen Kreis von — nach +, 

4 ee 0 p) so würde ein Tal verloren gehen, ohne daß H,_ı 
das Vorzeichen geändert hätte, |bewegen wir 

a de FTIR en dagegen zuerst den kleinen Kreis von + nach —, 

-_ a (RER — würde ein weiteres Tal im Inneren entstehen. 

“ Abb. 7. a) nicht mögliche einzahnige und b) mögliche Beides ist aber mit dem Satz von Orlando 


Raupen-Bewegung eines zweizahnigen Tales . Sue 
u rg ; nicht verträglich. 


Dagegen können die zweizahnigen Täler durch Verbreiterung und nachfolgender Zusam- 
menziehung, also durch die bereits oben erwähnte Raupenbewegung, beliebig um eine gerade 
Zahnzahl weiter rücken. Dies zeigt Bild 7b, worin zunächst links ein weiterer großer Kreis von 
+ nach —rückt, also das R-Tal um zwei Zähne verbreitert, und dann rechts der große Kreis von 
— nach + rückt. In diesem Falle wäre auch der Doppelschritt ein möglicher Grenzfall. Wir 
vermerken wieder, daß das Weiterwandern eines geraden Tales um gerade Zahnschritte nichts 
an der Zahl der Vorzeichenwechsel in der Hauptfolge oder Nebenfolge der H-Determinanten 
ändert, wohl aber wieder dazu führen kann, daß die Zeichenwechsel in beiden Folgen beliebig 
weit auseinander liegen. 

Wenn das Weiterwandern der beiden Täler amrechten Ende der H-Kette 4 positive Werte 
zurückgelassen hat, so kann ein drittes Tal entstehen, wie das in Bild 8 gezeigt ist. Damit ergibt 
sich als neue Möglichkeit der Zahnaustausch zwischen geradem und ungeradem Tal, wie er in 
Bild 9 eingezeichnet ist, und der nur dazu führen kann, daß das vorerst gerade Tal ungerade 
wird und das ungerade Tal gerade. Dieser Wechsel ändert aber nichts an der Zahl der Vor- 
zeichenwechsel in der Haupt- oder Nebenreihe der H-Determinanten. Dagegen zeigt Bild 9 
wieder deutlich, wie weit die beiden Vorzeichenwechsel in beiden Folgen voneinander abrücken 
können. 
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Die Entstehung weiterer Täler kann nichts Neues mehr bringen. Es können nur gerade 
oder ungerade Täler entstehen, deren Wanderungsgesetze wir bereits kennen: es können nur 
Zahnaustausche zwischen geraden Tälern oder — was einfach die Umkehrung dieses Vorganges 


+ 0-0 Oo OS G)EORE) o 


RER OL@)ECE@) euer 
= OÖ o BR Ze (6) SEE 70 
[») 
+ m —— + 
R NORA P AR WATER AI 
_ —x x ne - De 
Bild 8. Entstehung eines dritten R-Tales Bild 9. Zahnaustausch zwischen geradem und 
ungeradem Tal 
ist — Zahnaustausche zwischen ungeraden Tälern, oder schließlich Zahnaustausch zwischen 


einem geraden und einem ungeraden Tal entstehen. Da aber alle diese Teilvorgänge nichts an 
den Zeichenwechseln in der Haupt- oder Nebenfolge der H-Determinanten ändern, ist damit 
der behauptete Satz bewiesen. 


Eingegangen am 13. Juni 1952. 


Greensche Funktionen für elastische Platten 
Von F. Schultz-Grunow, Aachen 


Es werden gut konvergierende Entwicklungen für die zu AAu gehörenden Green schen Funktionen ange- 
geben bei einfach zusammenhängenden Bereichen beliebiger Umrißform und zwar für Randbedingungen, wie 
sie bei elastischen Platten auftreten. 

For simply connected regions of arbitrary shape, and for boundary conditions that occur in the theory of 
elastic plates, rapidly converging series eupansions of Gr ee n’s functions for AAu are given. 


Pour des regions simplement coherentse de forme arbritraire, sp&cialement & des conditions marginales, 
se presentant & des plaques &lastiques, des developpements bien convergents pour les fonctions de @reen appar- 
tenant a AAu sont indiques. 


Ana pyueunü Ipuna, orHocammxea x AAu, NAMOTCH XOPOIMO CXONAIINHECH PA3BIIOKCHUA MI 
ONHOCBASHEIX O6NACTeH IIPOHSBONBHOTO KOHTYPp& HA cAIyyah Kp&eBbIX YCIIOBHÜ, BO3HHKAIMIIMX 
y yupyrux IIIACTHHOK, 


Allgemeines 


Das Singularitätenverfahren für Rechteckplatten [1,2] versuchte bereits Greensche 
Funktionen herzuleiten, um der modernen Fragestellung nach den Einflußflächen, welches die 
technische Bezeichnung für Green sche Funktionen ist, entgegenzukommen. Es leidet aber 
an einer ungenügenden Erfüllung der Randbedingungen, weil die benutzten biharmoni- 
schen Funktionen längs des Randes beliebig oft ozzillieren, während dies bei den Funda- 
mentalintegralen nicht der Fall ist. Besondere Beachtung scheinen die Ecken der Berandung 
zu verdienen, wie der Verlauf der Spannungstrajektorien [3] zeigt. Hier werden dagegen die 
biharmonischen Funktionen, die dem Fundamentalintegralzu überlagern sind, aus harmonischen 
hergeleitet, wodurch funktionentheoretische Hilfsmittel anwendbar werden, die die exakte Er- 
füllung der Randbedingungen in den Ecken und ihrer Umgebung erlauben. 


Die zu AAu gehörende Greensche Funktion 


Die Greensche Funktion @ (&,, 0; & n) mit den laufenden Ortskoordinaten &, n und 
den Koordinaten &, 7, eines Aufpunkts in der Plattenebene setzt sich zusammen aus dem 8 - 
fachen Fundamentalintegral mit 


®?=(E —&)”+(m N VE (1) 
und den biharmonischen Funktionen F,, (&, 705 & N) 
G (EIR No> = N) — 0? In D-- = 18 (Eos No5 &, n) NETTE TEN (2), 


wo SF, durch die Forderung bestimmt ist, daß @ gewissen homogenen Randbedingungen 


m 
zu genügen hat. Als solche kommen in Frage, wenn w die Durchbiegung, n die Normale, 
ti die Tangente des Randes, u die Querkontraktionszahl bedeuten, 
ou 


an Eh N A 63) 
w—=0, Sri (3) 
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bei eingespanntem Rand, 


0°u u 
u=0% mil RE 
bei Navierscher Lagerung und 
ou ou u ou 
—— _——_ ar 2 U) — EN Eee 
Oma Be De na 6) 
am freien Rand. Z. B. muß im Fall der Randbedingungen (3) F den Randbedingungen 
0 0 
PTR . 2 Me (92 Kae ENG 
DIR eng; San 37 (0° In 0) (6) 


M m 


genügen. 
Bei homogenen Randbedingungen läßt sich die Lösung von 


Adu=p(E,n) 


stets darstellen in der Form 


Wir 
u (&0, 70) = | | Glenn) PEN) dEdn Me et) 
r 


wo p (&, n) stetige Ableitungen 1. Ordnung besitzen muß. 

Im Fall der Platte mit der Dicke A und mit dem Elastizitätsmodul E bedeutet u die mit 
dem Faktor Eh?/12 (1 — u?) multiplizierte Durchbiegung, 9 (&, n) die auf die Flächeneinheit 
bezogene Flächenbelastung der Platte im Punkte £,n und die Singularität des Fundamental- 
integrals im Punkte &,n, bedeutet eine dort wirkende Einheitskraft Q = 1. Um letzteres zu 
zeigen, denkt man sich um den Kraftangriffspunkt als Mittelpunkt einen beliebig kleinen 
Kreisschnitt des Halbmessers a gelegt, damit das Fundamentalintegral in jeder gewünschten 
Größe die Funktion F überwiegt. Für die im Schnitt wirkende Querkraft q je Längeneinheit 


des Schnittes gilt g= — Au. 


Die totale, im Schnitt wirkende Querkraft ist, wenn der Umfangswinkel mit bezeichnet wird, 
2n 27 


l 0 
fesar = | 2 (a0) up. 
() ö 


Hier für w das Fundamentalintegral eingeführt, erhält das Integral auf der rechten Seite den 
Wert Eins. Das auf der linken Seite stehende Integral ist aber aus Gleichgewichtsgründen der 
wirkenden Einzelkraft © gleich, so daß man tatsächlich@ = 1 erhält. @ (&,, no; & n) bedeutet 
also die Biegefläche unter dem Einfluß einer im Punkte &,, 7, wirkenden Einheitskraft. Die 
technische Bezeichnung für diese Fläche ist Einflußfläche der Durchbiegung wegen der Vertausch- 
barkeit des Kraftangriffspunktes (£&,, 70) mit der Stelle &, n, an der die Durchbiegung bestimmt 
werden soll, infolge der Symmetrie: @ (&,765&n) =@(&,n; &, 0) oder physikalisch aus- 
gedrückt, als Folge des Satzes von Betti. Dieser Satz über die Gleichheit der Ergänzungsarbei- 
@ (&0» 005 8 N) 
8% 


erzeugte Durchbiegung bedeutet, die von der Flächenlast 9 (&, n) erzeugte Durchbiegung % 
1 
U (&> 70) | [e% No5 5) pP (6, m) dE dm, 


welcher Ausdruck mit (7) übereinstimmt. 


Sobald @ für einen Aufpunkt ermittelt ist, so ist damit die Durchbiegung in ihm für 
irgendwelche Lastfälle bestehend aus Einzellasten oder der obigen Einschränkung genügenden 
Flächenlasten gegeben, denn für diese läßt sich die Integration in (7) elementar durchführen 
und bei Einzellasten kann man stets auf die im Aufpunkt wirkende Einzelkraft zurückführen, 
wegen der erwähnten Vertauschbarkeit der Kraftangriffsstelle mit der Stelle, an der die 
Durchbiegung angegeben werden soll. Ist nicht die Einflußfläche der Durchbiegung, sondern 
etwa die des Biegemoments gesucht, so kann auch diese unmittelbar hergeleitet werden, denn 
in (7)ist nach den Aufpunktkoordinaten zu differenzieren, was unter dem Integrationszeichen 
geschehen kann, so daß auch Greensche Funktionen für die Schnittreaktionen bestehen. 
Die Randbedingungen sind für diese, als höhere Ableitungen von u die gleichen wie für « [1]. 


ten liefert, wenn die von der im Punkte &,n. wirkenden Einheitskraft @& 


G 


Pe 
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Konstruktion der Bipotentialfunktionen 


Die dem Fundamentalintegral zu überlagernden Bipotentialfunktionen F,, werden in 

der Form 
Il OH) an (8) 

aus Potentialfunktionen gebildet, die in der Kreisebene (z-Eebene), auf die sich die Platten- 
ebene (S-Ebene) stets konform abbilden läßt, bestimmt werden. o ist der Fahrstrahl in der 
&-Ebene, wie er durch (1) gegeben ist. Diese Darstellung biharmonischer Funktionen ist 
erlaubt, wenn der Bereich den Nullpunkt enthält und der von ihm ausgehende Fahrstrahl o 
den Rand nur einmal schneidet. Die Funktionen y könnten in die Plattenebene transformiert 
und dort dann die Randbedingungen erfüllt werden. Jedoch ist es bei den i. A. recht ver- 
wickelten Abbildungsfunktionen zweckmäßiger, die Transformation nicht erst auszuführen, 
sondern umgekehrt die Randwerte in die Kreisebene so zu übertragen, daß zugeordnete Rand- 
punkte gleiche Werte besitzen, und die Randbedingungen in der Kreisebene zu erfüllen. 

Um die Randwerte der Ableitungen zu befriedigen, müssen die Ableitungen der Bipoten- 
tialfunktionen F,, selbstverständlich in der &-Ebene gebildet werden, z.B. die Normalab- 
leitung 


OF, , dm (%, %) vu „2 O%m (&, %) N 00° C 
Ey RS ISEREER en 9 2 er): 

on Um on = n 0 on r 05 Ym. (&, Yy) on ( ) 
wo n die Randnormale der ö-Ebene bedeutet. Die Ableitungen der Funktionen y,, (2, y) in 
der &-Ebene lassen sich mit der Beziehung 


dw„(2) _ dw, (z) de 
Bed de 


aus den Ableitungen in der z-Ebene errechnen, wo w die komplexe Funktion bedeutet, deren 
Real- oder Imaginärteil % ist. 

Von besonderem Interesse für die Plattenbiegung sind polygonale Berandungen. Hier 
bringt die konforme Abbildung auf den Kreis den Vorteil, daß die Unstetigkeit der Richtung 
in der Ecke verschwindet. Die Ecke wird zu einem Häufungspunkt äquidistanter Randpunkte. 
Für die Abbildung des Plattenbereiches auf das Innere des Einheitskreises besteht die 
Schwarz-Christoffelsche Formel 

Ne in deutet), 
dz K 
wo A eine Konstante und 2x die komplexen Koordinaten der Randpunkle des Kreises sind, 
deren Bildpunkte in der &-Ebene Eckpunkte der Berandung sind. 9x ist eine positive Zahl und 
zwar ist, wenn der am Kreis bestehende Randwinkel x bei der Abbildung um das g-fache ver- 
kleinert wird, so daß der Eckenwinkel gr beträgt, 


ri) 


geiz. 
1 : : 
Z. B. ist an einer rechtwinkligen Ecke g -. und daher 9x = 5 ; bei einem Eckenwinkel von 
Pa istg ; En ? = 
eg ET 


Um die Art der Singularität des Randwertverlaufes des Fundamentalintegrals in der 
Kreisebene zu erkennen, sei die Abbildung in der Nähe des Bildpunktes 2x einer Ecke betrach- 
tet. Hier läßt sich (11) schreiben 

de 

a 
wo V(z) im Punkt 2x nicht Null und nicht Unendlich wird, so daß sich V(z) in eine Potenzreihe 
um den Punkt 2x entwickeln läßt. Hiervon genügt für die Untersuchung der Singularität das 
konstante Glied a, so daß sich 


Gen Ev). 


ne 5 


dz 


ergibt. Der Kürze halber sei 2% = 1 gesetzt, d.h. der Bildpunkt des Eckpunktes liege auf der 
reellen Achse, was keine Einschränkung bedeutet. Dann gilt für Randpunkte, wenn a der 
Winkel des Pfeiles z mit der positiven reellen Achse ist, 

& 


—g ’ Ike z 
(2m?) en se rn (13), 
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Mit 
= = = (2 sin 2] e- „le 
ergibt sich 

2 


wo der Absolutbetrag von d£ das Bogenstück ds auf dem Polygonrand bedeutet. Für die Um- 
gebung der Ecke besteht also, wenn s von der Ecke aus gezählt wird, 


2|- &\9E 
a -@l2sng) eo. (ana: 


woraus man die Zuordnung der Randpunkte der &- und 2-Ebene in der Nachbarschaft der Ecke 


un) 19x 
68, — Mt ae re a 
gr 149) 


erhält. 


I 
EWR 


Bild 1. Zuordnung der Randpunkte 


Um den Randwert verlauf des Fundamentalintegrals in der Kreisebene in Eckennähe zu 
erhalten, liest man in Bild 1 @ f 


ee LE . (16) 


e=b! —2cs +2=b? —2c 
I IE (eg): 


ab. Hiermit ergibt sich 


w|- 


2 
eelno= BZ: b°—2c(l+In (in ve] 


und mit (15) 


1 17 a: 
2] nn. b2 2 2 0 IK 0 2(1—g7) 
0 ne=;| In b 2c(l land), — Pi rim ale zur Bl Al2) 


Ferner besteht ein Randwertverlauf einer Ableitung, z. B. bei eingespanntem Rand der Nor- 
malableitung. Es ist (s. Bild 1) 


00? 
on — 2 d . . * * * . ” * * ” . . * * * * (18). 
Hiermit 
do? In o R a Er 
Tr = Hm) ram ee a Ba LASER: 


wobei wiederum Glieder einschließlich 2 (1 — gx)-ter Ordnung berücksichtigt sind. 
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Mit (17, 18) sind die Randsingularitäten der Fundamentalintegrals bekannt und es lassen 
sich nun sofort Funktionen angeben, die gleiche Singularitäten aufweisen und deshalb geeignet 
sind, die Randbedingungen in der Ecke K exakt zu erfüllen. Dazu werde die Funktion 


v=y+ip—lar —2)HE .....2.2.2.2.0. . (20a) 


betrachtet. Es gilt am Kreisrand, wenn wieder 2x = 1 eingeführt wird, beispielsweise für den 


Realteil 
BRR\ SIR) 3 & 
v—(26n3) cos ar +3)« |. 


Die ersten beiden Glieder der Reihenentwicklung dieses Ausdrucks nach a ergeben 


= 3 1— 3 
v=a' "Rs. m(l —g)—a E sinn (l — gg) 2.0.5 (206). 


Ferner interessiert eine Ableitung dieser Funktionen am Rande und zwar bei der als Beispiel 
gewählten Randbedingung der Einspannung die Normalableitung. Wegen 


BE. ...........@ 
on 92 orlde 
ergibt sich mit (15) 
Buell 0E% Ic DOM 
ne % cos | 2) #30 n 
Hier ergibt die Reihenentwicklung nach « 
GRAUEN, 3 IE\ .; Be 
re rn BE 75981 (i , sinn ge (22) 


Führt man (16, 18, 20b, 22) in (8, 9)ein, so erkennt man, daß die aus y mit geeigneten Werten 9x 
gebildeten Bipotentialfunktionen in der Ecke die gleichen Singularitäten in den Randwert- 
verläufen besitzen wie das Fundamentalintegral, so daß mit diesen Bipotentialfunktionen die 
Randbedingungen in der Ecke und deren Umgebung exakt erfüllt werden können, d. h. daß die 
Glieder mit a-Potenzen <1 in den Reihenentwicklungen des Fundamentalintegrals und dieser 
überlagerten Bipotentialfunktionen zu Null gemacht werden können. 

Zur Erfüllung der Randbedingungen an weiteren Randpunkten eignen sich Bipotential- 
funktionen, die sich in der Form (8) aus komplexen Funktionen mit Nullpunkt im Kreisinnern 
herleiten lassen und die keine Unendlichkeitsstellen im Kreisinneren und auf dem Rand haben, 
also etwa 


Die aus solchen komplexen Funktionen hergeleiteten Bipotentialfunktionen besitzen gleiche 
Randsingularitäten wie das Fundamentalintegral, was bei Funktionen der Form 


Fn = Cm 0” cos me + Om or+?2cos me 
nicht der Fall ist. Außerdem läßt sich 
F=%+0(0. Wa er 1 Pr . (24) 


verwenden. 


Die quadratische Platte 


Die Platte sei beispielsweise als eingespannt vorausgesetzt und es wirke eine Einheits- 
kraft in der Plattenmitte. Die Schwarz-Christoffelsche Formel (11) lautet hier mit 


1 
IK = 5 
d + 
u 
dz/a a 
wenn die Eckpunkte des Quadrats (£-Ebene) auf den Koordinatenachsen und die Mittelpunkte 
von Quadrat und Kreis (2-Ebene) in den Koordinatennullpunkten liegen. c bedeutet die halbe 


Seitenlänge des Quadrats, «a den Kreisradius, x ist ein noch zu ermittelnder Zahlenwert. Der 
Kürze halber wird an Stelle von ö/c und z/a stets £, z geschrieben und zwar werden sich die 
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Koordinaten im folgenden ausschließlich auf Randpunkte beziehen, weil es sich hier insbesondere 
um die Erfüllung der Randbedingungen handelt. Dann gilt 

1 ..n—4& 
dd TE 5) 


7, + (@sin 20) e ee a 2: 


und 


ug 
dE dEde ; ur 

—_ — sin 2 Be . 
A PErE x (2 sin 26) e 


Ist s/c bzw. s die von der auf der positiven reellen Achse liegenden Ecke aus gemessene Längen- 
koordinate auf dem Rand des Quadrates, so gilt 


ds de -4 
Sn ar 2 
75 12 x*(2sin 20) 
und 
ya 10 Bram a (26). 
Y2sin 2 


Dem Wert s = 1 entspricht x = n/4. Da das Integral zwischen den Grenzen 0 und n/4 den 
Wert 1,854074/2 hat, so ergibt sich x = 1,07875. Für die Nachbarschaft der Ecke liefert (26) 
folgende Zuordnung der Randpunkte 


se Hya.. one ee 

In der Nachbarschaft der Ecke ergibt sich mit (27) (Vgl. Bild 1) 
elite 2— ante Mn ......(2) 
ing =m2 —x(l +1n2)ar +% (2 HE jn OD). 2 0 02 SEE a 23 
elng=1+m@®=1+4m2—xai. abet un Aw eaeheind Hau 


Hiermit sind die Randsingularitäten des Fundamentalintegrals in der Kreisebene bekannt. 
Zum Ausgleich der auf der positiven reellen Achse liegenden Randsingularität wird 


gemäß (20a) mit gx =, die Funktion 


wıı — ne) 


gebraucht und zwar wegen der vorliegenden Symmetriebedingungen ihr Realteil 


1/2 
ln) lt): 


im Bereich 0 <a <2x. Ihre Normalableitung am Rande ergibt sich mit (21, 25) als Realteil 
der Normalableitung von w,, zu 


Op 1 (sin2%\ 142 3 7 : 
FR Eee (1° u ie ee. (228 
sin > 


Ferner wird die Funktion 
0, = (1 2) 


benutzt, um auch in der Nachbarschaft der Ecke die Randbedingungen exakt zu erfüllen, 
und zwar wegen der vorliegenden Symmetriebedingungen ihr Realteil 


nA LE ) 3 y 
%ı =2 (ir 2) cos (7° dr J) ee Se 
genommen. Ihre Normalableitung am Rande lautet 
a RR. Hall 28 5 
7 -— [sing sin 2 cos(4 +4) eier reee 


ebenfalls für den Bereich 0 <a <2n. 


G 


o 
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Die im Randpunkt & = 0 bestehenden Reihenentwicklungen lauten 


1 er & oyıı 1 3 ) 
= u ed 1 € =. ; 
Yı N #7); an Sr 1 7); 
Y%ı=0; Opa 2 


on \2% & 
einschließlich der Glieder in a. 


Die weiteren Potentialfunktionen y,x, Ysx, die für die übrigen Ecken gebraucht werden, 
lassen sich leicht aus (31, 32, 33, 34) durch Drehen des Koordinatensystems erhalten, wobei 
allerdings der Geltungsbereich so zu wählen ist, daß man im gleichen Riemann schen Blatt 
bleibt. Z. B. ist der Geltungsbereich der Funktionen %5; % für die auf der positiven imagi- 


nären Achse liegende Ecke: = Zn <on, so daß die Reihenentwicklung für den auf der posi- 
tiven reellen Achse liegenden Randpunkt nach &« — 360 und nicht nach & vorzunehmen ist. 
Das gleiche gilt für alle übrigen Funktionen. Auch für %,,, % Kann wegen ihrer Symmetrie zur 
reellen Achse « durch 360 — a ersetzt werden. Man hat die folgenden Ausdrücke, neben welchen 
der Geltungsbereich angemerkt ist 


NUR ER ZEN LÄLS Eur (0061 Ir 9 7 5) 
m 205 7)| (dt) Dr 
a AU 6 
ae NT REN nd 1 
Yı3 | (9 3] sin 1 am 5% 


Ya —|2sin u, 
2 4 


mit den Normalableitungen 


a — 1/2 5) = 
nn — m 2 a)? si © — al sin (f IR u ) 
En — 1/2 3 
is je : (sin 2 a)t/2 in [E = 2) an (e Be ) 
Oyıa 1 A. 3 2. 13 a 
Mm 2x 5 Di BETT gr me . 
on 0) „(sin Z &) /2 | sın 5 1 7L sin n n ; r 


Die Reihenentwicklungen für die auf der positiven reellen Achse liegende Ecke lauten, ein- 
schließlich der Glieder mit der ersten Potenz des Winkels, wenn 360 —ıa — ß gesetzt wird, 


4 IT 1 TT = IT ya 1 | E 2) By 
—= —Y2jcs——+-—[cos— + sin i ; =—; sin — | ß!/? 
Bi | 8 A| B ;) on TEE: 


ae. Oy il 
Ya —V2; oz a BR (35). 
2 1 4 
Yu=— y2 E = — ze = sin a ß s ei — — (Mn :) ß}/® 
C y2 x 
Ferner hat man 
Oi & 7.\ ]3/2 3 u 
Ya sın ZEW/ U cos BUT ’ 
ge rer Snap ınlan|[% z\ |"? . 5 ar 
en (en 26) sin|>—7 cos 73 
3/2 
Ya = | 2sin ei cos men : 
2 2 4 2 
1/2 5 IT 
Yan ze (sin 2 0)" I 2 _ 3)! cos| 7 a — A 
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A age JE FROM e A R 
DR ae er 
Mal ar 9,8 ml 1% 3 nr 5) 2 
Ger == = (sin 2 «) sın I A cos A & BE 
und die Reihenentwicklungen 
3 7 u 0 3 m 
Bose in ; Va) Mi ae Ro 
We 2 I a] (eo 8 + sin ;)?| ; Mr 2 21/% sin R ß 
0 3 
N — 23/2; ee en ER 91/2 pr? (36). 
a ae 3 7 Ei, Ä U ae 
Pig I R 7 (es F sin s)ßl; nl sin R DW 


Man erkennt, daß sämtliche Funktionen yı x, Yax geeignet sind, die Randsingularitäten 
des Fundamentalintegrals, die durch (29, 30) gegeben sind, auszugleichen. Bildet man jetzt mit 
den Reihenausdrücken die Bipotentialfunktionen F,x, F5x gemäß (8) mit je gleichen Kon- 
stanten (’, C’’ wegen der bestehenden Symmetriebedingungen, sowie die Normalableitungen 
nach der Anweisung (9), so kann man zunächst in den Ecken und ihren Umgebungen die Singu- 
laritäten eliminieren, indem man unter Berücksichtigung von (29, 30) 


In2 —x(1 + In 2) 81? ee +1m2)P+ Dr + Fr) =0 
K=1 


4 

1 +4In 2 — RD + > (Er + F,x) —=( 

K=1 
setzt und die Konstanten 
0’ C’’ so bestimmt, daß die 
Glieder mit gleichen Expo- 
nenten von ß je für sich Null 
werden. So ergeben sich vier 
lineare Gleichungen zur Er- 
mittlung der Konstanten 
Ci, Ci, 03; 0%. Man erhält 
die Werte 


; Ä — er C4= 0,8187 
En IV(100- fach) 2 2 
en r — 0,3534, 


I 
C3—= — 1,2059, 
2= 0,5247. 


Den Verlauf der Randwerte, 

# die nach Überlagern der 

N 2 TI are unkioun F,x, Fax über 
DB 1 02 0 Won une s das Fundamentalintegral 
übrig bleiben, zeigen bis auf 


eine Konstante die Kurven Il 


IT (10- fach) 


Bild 2. Randwertverläufe der Durchbiegung bei verschiedenen Näherungen 
I Randwertverlauf des Fundamentalintegrals 


II Rand wertverlauf nach Ausgleich der Randsingularitäten in Bild 2; 3: Die Kurven I 
III Zusätzlich zu II Erfüllung der Randbedingungen beix = 22,5und 45° . B ’ 
IV Zusätzlich zu III Erfüllung der Randbedingungen bei& = 10° zeigen zum Vergleich die 


Randwertverläufe des Fun- 
damentalintegrals. Es bestätigt sich, daß die Randbedingungen in der Ecke und in ihrer 
Nachbarschaft exakt erfüllt werden können, denn die Randwerte sind dort Null und sie be- 
sitzen eine horizontale Tangente. 

Um die restlichen Randwerte zum Verschwinden zu bringen, werden gemäß (8, 23) Bi- 
potentialfunktionen der Form 
Fn= Cm rt"? cos ma + On o?rm cos ma 
überlagert. Wegen der bestehenden Symmetrien kommen für m die Zahlenwerte 
mA, 8, 12,23> 
in Frage. Ferner wird die durch (24) gegebene Funktion F, überlagert. 
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Die Normalableitungen der Potentialfunktionen r” cos m a am Rand, aus denen die Funk- 
tionen F,, gebildet sind, ergeben sich auf dem oben beschriebenen Weg zu 
0 m 
—ym ner et 
An EM ( 
und es bestehen die Reihenentwicklungen um ß =0 (ß = 360 —.a) einschließlich der in ß 
linearen Glieder 


2 sin 2 )U/2 cos ma 


1) Mm 
om BL; En cos mi = — 2 pi", 

Zunächst wurden zwei Funktionen F,, mit m = 4,8 zu den Funktionen F, Fıx; Fax hin- 
zugenommen und damit in den den Winkeln & = 0, 22,5°, 45° entsprechenden Randpunkten 
des Quadrats die Randbe- ar 
dingungen erfüllt. In jedem On 
dieser drei Randpunkte be- 
stehen zwei Randbedingun- 
gen, außerdem in der Ecke 
(«=0) die vier Bedingungen, 
daß die mit 1/2 und ß behaf- 
teten Glieder des Randwert- 
verlaufes und seiner Normal- 
ableitung verschwinden. Die- 
sen 10 Bedingungen stehen 
10 freie Konstante Ü gegen- 
über. Dabei zeigt sich, wie 
sich auch allgemein nach- 02 
weisen läßt, daß zwei für die 
Ecke bestehendeBedingungs- -94 
gleichungen linear abhängig 
sind, nämlich die Bedingung 
des Verschwindens der mit 
ß!/? behafteten Glieder des 
Randwertverlaufes und des #8 R 
konstanten Gliedes der Nor- Bild 3. Randwertverläufe ae beiden gleichen 
malableitung. Es bleibt also 
eine Konstante frei. Sie wird so bestimmt, daß die Normalableitung möglichst mit gleichen 
Amplituden um den Randwert Null oszilliert. Es ergeben sich so Randwertverläufe, die die 
Kurven Illin Bild 2, wiedergeben. 

Endlich wurde noch die nächste Funktion F,, mit m = 12 hinzugenommen, womit zu- 
sätzlich die Randbedingungen auch bei & = 10° erfüllt wurden. Die Zahlenwerte der Konstan- 
ten sind in diesem Fall 


FE 
150 20° 22,5° 0° 35° ug° 4g0Z[ia5 @ 
u50 —n 

02 


II (10 fach) 


IF (10 fach) 


0,6 


01 = 0,633703 , 1= — 0,445908, (&= — 0,727508, (= 0,409825 , 
O4= — 0,067793,  Ci= 0,021776, C3= 0,002753 , 03 = — 0,000908 , 
O2 = 0,007392 , Cie = — 0,00659 . 


Die bei dieser Überlagerung noch verbleibenden Randwerte zeigen die Kurven IV in Bild 2, 3. 
Der Vergleich der Kurven II, III, IV mit den Randwert-Kurven I des Fundamentalintegrals 
ergeben, daß diese Berechnungsmethode tatsächlich gute Konvergenz besitzt. Der Maximal- 
wert des verbleibenden Randwertes beträgt bei der genauesten Näherung beim Randwert selbst 
0,2%, des Maximalwertes des Fundamentalintegrals und 0,4% bei der Normalableitung. Ohne 
die Funktion F,, beträgt die Genauigkeit 1% und 4,9%, welche i. A. bereits ausreichen dürfte. 

Nachdem die einzelnen Funktionen ermittelt sind, läßt sich die Gr ee nsche Funktion @ 
für verschiedene Punkte der Kreisebene berechnen und dann in die Bildpunkte der Platten- 
ebene übertragen, wobei die von Sch warz [4] durchgeführte Berechunng der Bildpunkte, 
die auf vier Stellen genau ist, verwendet werden konnte. 

Die Biegefläche, die man schließlich so für eine in Plattenmitte angreifende Einzellast 
Q =1 erhält ist in Bild 4 für die beste Näherung dargestellt. Sie stellt die Green sche Funk- 
tion @ für die Durchbiegung dar. Ihre Ordinatenwerte sind in Verbindung mit Bild 5 durch 
Tabelle I gegeben. 

Um;die Genauigkeit dieses Verfahrens zu prüfen, wurde ein Vergleich mit den genauesten 
Werten durchgeführt, die für die Durchbiegung bekannt sind. Allerdings beziehen sich diese 
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Werte auf eine gleichmäßig verteilte Last. Dieser Fall konnte aus dem hier behandelten leicht 
hergeleitet werden durch graphische Auswertung des Integrals in (7) mit p = konst. Nach der 
eingangs gegebenen Definition der Größe «gilt, wenn f die Plattendurchbiegung ist, 

u=7,B 
mil 


B Eh 


RAS 


® 


Bild 4. Greensche Funktion für die Durchbiegung für Belastung durch Einzelkraft. Aufpunkt in Plattenmitte 


Es wurde, wenn a die Seitenlänge der Platte ist, 

u= p a* 0,001266 
erhalten. Genau den gleichen Wert errechnete Tölke [5], während Nadai[6] den Zahlen- 
wert 0,001268 angab. Es besteht also vollständige Übereinstimmung mit den bekannten ge- 
nauen Werten. 


Diagonale 


Bild 5. Zur Zuordnung der Tabellenwerte 


TabelleI. Ordinatender Einflußflächen 


0—0 0,563714 1—1 0,475004 2—2 0,326644 3—3 0,174829 
0—1 0,513594 1—2 0,393374 2—3 0,241445 3—4 0,103881 
0—2 0,418374 1—3 0,292384 2—4 0,155998 3—5 0,047008 
0—3 0,311264 1—4 0,193711 2—5 0,083951 3— 6 0,011347 
0—4 0,209216 1—5 0,110912 2—6 0,032928 3—7 — 0,000673 
0—5 0,123085 1—6 0,051171 2—7 0,005339 

0—6 0,060168 1—7 0,018983 2—8 0,000232 4—4 0,048627 
0—7 0,021998 1—8 0,007542 4—5 0,014422 
0—8 0,004933 1-9 0,000196 4—6 0,001506 
0—9 0,000154 j 


0-10 0.000037 5—5 0,000097 
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Es sei auch der für den ursprünglich behandelten Lastfall, Einzelkraft Q in Plattenmitte, 
erhaltene Wert der Durchbiegung in Plattenmitte besonders angeführt, da hierfür auch noch 
keine genaue Berechnung vorliegt. Der Tabelle Ientnimmt man 


u =Q a? 0,00561 . 


Die numerischen Resultate verdanke ich meinem Mitarbeiter, Herrn P. Schmitz. 
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Eingegangen am 1. Juli 1952. 


Elasto=plastische Reflexion eines Stabes *) 
Von Hans Richter in Haltingen/Baden. 


Es wird der elasto-plastische Reflexiensvorgang und der nach der Reflexion verbleibende Schwingungs- 
zustand für einen Stab endlicher Länge untersucht, der mit genügend hoher Geschwindigkeit gegen eine feste 
Wand stößt und dessen Material keine Verfestigung, jedoch Relawation besitzt. Die Abhängigkeiten des 
Restitutionskoeffizienten und der Reflewionsdauer von Stablänge und Geschwindigkeit werden berechnet. 


A bar of finite length is thrust against a fived wall with sufficiently high velocity; the material shows no 
hardening but relawation. T'he author investigates Ihe elasto-plastic phenomenon of reflecton and the state of 
oscillation remaining after thereflection. It is shown, how the coefficient of restitution and the time of duration 
of the reflection depend on the length of the bar and the velocity. 


Une barre de longueur finie, le material de laquelle ne possede pas d’&crouissage, mais se reläche, pousse ü 
vitesse assez haute contre un mur dur. L’auteur ewamine le proced£ de r&flexion elasto-plastique ei l’&tat d’oscil- 
lation restant apres la reflexion. Les dependances du co£fficient de restitution et de la durde de reflexion de la 
longueur de la barre et de la vitesse sont calculees. 

Crep:keHb KOHeYHOH MJIUHHBI YAAPAETCH C NOCTATOYHO BEICOKOH CKOPOCHIO O IKECTKYIO CTEHKY; 
MATEPHAJI CTEPIKHA He UMEET 3ATBEPAECHHA, ODIAAAET, ONHAKO, persakcammei. B cTarbe nechne- 
AyOTcA yupyre-macTu4yHbIa IPONMecc OTPAsKEHHA HM KOJEÖbATENIBHOE COCTOAHHE, OCTAWINEECH 
TOCJIO OTPAKECHHA. BeryucasioTtcH BaBHCHMOCTH koshphunuuenHrta BOCCTAHOBJIEHHA H AIIUTEJIBHOCTH 
OTPA3KeHUA OT UIHHBI CTEPIKHA U CKOPOCTH. 


1. Problemstellung 


Ein spannungsfreier Stab konstanten Querschnittes bewege sich mit konstanter Ge- 
schwindigkeit v, in negativer &-Richtung auf die feste Wand & = 0, die er im Zeitpunkt i = 0 
trifft. v, sei so groß, daß die Fließgrenze o, beim Aufschlag überschritten wird. Die Aufgabe, 
den anschließenden Reflexionsvorgang rechnerisch zu verfolgen, wurde bereits 1949 durch 
V.S. Lenskijt!) unter der nicht unanfechtbaren Voraussetzung behandelt, daß das Stab- 
material auch für dynamische Vorgänge dem statischen Prandtlschen Spannungs-Deh- 
nungs-Gesetz gehorcht. Er erhielt hierbei hinter der vorauseilenden elastischen Stoßwelle 
konstanter Intensität noch eine plastische Welle von geringerer Fortpflanzungsgeschwindig- 
keit. Aus der akustischen Durchrechnung der zahlreichen Frontdurchkreuzungen konnte 
Lenskij Aussagen über die Reflexionsdauer, sowie über die Ausdehnung und Stärke der 
notwendig stückweise homogenen Deformationsgebiete gewinnen. Hier soll nun die gleiche 
Aufgabe für verfestigungsfreies Material behandelt werden, das bei dynamischer Beanspru- 
chung Relaxation zeigt, so daß unmittelbar hinter dem elastischen Wellenkopf plastische 
Defor mationen auftreten, die die Stoßwellenintensität im Laufe der Zeit schwächen. Wie bei 
Lenskij geschieht die Rechnung in akustischer Näherung. 


2. Die Grundgleichungen 
Hat ein Stabelement, das sich bei { = 0 an der Stelle & befindet, zur Zeit {> 0 die Span- 
nung o(z, t) und die Verschiebung «(z, t), so lautet die Bewegungsgleichung 


our =): . . . . . . . . . . . . . . (2.1), 


*) Auszugsweise vorgetragen auf der Jahrestagung 1952 der GaMM in Braunschweig. 
1) V,S.Lenskij, Über den elasto-plastischen Stoß eines Stabes auf ein starres Hindernis. Priklad. 
Mat. Mech., Moskau, 13 (1949), S. 165— 170. 
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wenn o die Dichte ist. Die erreichte totale Dehnung e = u, ist bei infinitesimalen Verzerrungen 


die Summe aus elastischer Dehnung &7, = 5 und plastischer Dehnung 


[07 
Ep = Un age 5 nn sage ae el 


wobei E der Elastizitätsmodul ist. &,; ändert sich zeitlich nicht, solange |o|<o, ist; für 
|o|> o, dagegen trete Fließen gemäß 
d j 
Hm = 0 — 0, sign o 
ein. Allgemein ist also 
ul %) —.sign a’ !maxX (0, of 0,) 7 m Fear Ve 2) 


Bewegungsgleichung (2.1)und Stoffgesetz (2.3) liefern das zu Grunde zu legende hyperbolische 
System partieller Differentialgleichungen für Spannung o und Geschwindigkeit v = u: 


0, —0,=0, | 
Eu — = sion co: max (0, |o| — o,) nn (A) 
Die zugehörigen Charakteristiken sind 
=, ett s)—const und n= 4 (et — a) = const bei c= n 
mit den Charakteristikenbeziehungen für stetige Wellen 
d(o + ocv) +2 sign o - max (0, |o| — o,) - al; — (für | = Bonstä a SAW): 


Wegen des Nichtschneidens von Charakteristiken derselben Schar pflanzen sich auch endliche 
Störungen längs der Charakteristiken fort und genügen den aus (W) folgenden Stoßfrontbe- 
ziehungen 


A(ofocv)=0 für Fronten „= const uhr age ehe: (DE), 


Die Wellenköpfe sind also stets rein elastischer Natur entsprechend der Tatsache, daß bei un- 
stetigen Änderungen keine Zeit für plastische Deformationen zur Verfügung steht, die erst im 
anschließenden Wellenzug einsetzen. (W) zeigt weiter, daß im elasto-plastischen Gebiete stets 
rechts- und linkslaufende Wellen gekoppelt auftreten. 

Im elastischen Gebiete |o|< o, führt (G) zur bekannten D’Ale mbertschen Lösung. 
Im elasto-plastischen Gebiete |o|> o, jedoch wird (G) durch 


EN 
ee 7. (1%) 
k Sr 


ei 
0 —- 0, SIgN Fr; RE BER (Ke+%n— 2%) 


gelöst, wobei x(&, n) bei der Normierung x (0, 0) = 0 als Lösung von %:, = x die quellenmäßige 
Darstellung mit Hilfe der Bess el funktion nullter Ordnung gestattet: 


€ Ast: n sich: 
ir } Jol2 iyyn)f Ey) Buch Jl2iyye)g’(n—y)dy mit f(0)=g(0)=0 . . (L.). 
Dabei ist (8, 0) = f(&) und 4(0, n) = g(n). (L,) in (L,) eingesetzt liefert Darstellungen für v und 
o, auf deren Wiedergabe hier verzichtet werde. Die spezielle Gestalt der Quellstärken f und g 


ist aus den Rand- und Anfangsbedingungen zu ermitteln. 


3. Der Stoß des unendlich langen Stabes 


Trifft der Stab zur Zeit £=0 mit der Geschwindigkeit v—= —v, auf die feste Wand 
& = 0, so entsteht ein Wellenkopf längs n = 0, der im & - t-Diagramm von Bild 1 vom Gebiete I 
mit o=0,v= — ins Gebiet n>0 führt. Gemäß (St) ist 


o&,+0)+ecW&5+N)+ul=0....... u EN. 
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Aus der Randbedingung v(0, + 0) = 0 folgt sofort 
(0, +0)= —opcvw. 
Damit überhaupt plastische Deformationen auftreten, muß also 
} o 
V>d, bei v,= 5 


sein, was im folgenden stets als erfüllt angenommen ist. Die dann 
längs n=-+0 gültige W-Beziehung führt unter Beachtung von 
sign o—=— 1zu 


d d 
TE o(£&, +0) 007506 +0)+2-[o(& #0)+0,]=0, 


was zusammen mit (3.1) liefert: 


EN) __ TRADE AR PR NEN 


V— ds ec — %,) 


Bei Übergang zu dimensionsbehafteten Größen bedeutet dies, daß 
die Intensität des Wellenkopfes vom Anfangswerte ocv, mit oc/2u 
als Dämpfungsfaktor bezüglich z exponentiell abklingt bis zur 
Fließgrenze o,. Dementsprechend kann die Materialgeschwindig- 
keit nur bei &=0 annulliert werden, während für 2>0 eine auf 
die Wand gerichtete Restgeschwindigkeit verbleibt, die bei 2 — 
bis zu %, — v, wächst. 


Bild 1 


Um nun die für n>0 maßgeblichen Quellstärken f und g zu 
erhalten, haben wir die Randbedingungen v(£=n) —=0 und (3.3) in (L,,) einzuführen und 
erhalten nach einfacher Rechnung: 


& — 

(0) Hai) + SA Ey) EN — 2 N] dy=0 bei z=f—9, 
& 

70-04 [Pro —ro1ae = N ae), 


& 
HaUERAUE: j LE) —2S (E) + FO] dE = — erde Ew1 


Dieses Gleichungssystem hat die eindeutige Lösung f(y) = g(y) = — un 2uye’, so daß 


wir nunmehr nach einfacher Rechnung erhalten: 


n 


re 1 eilEm) 20 | @iyEN er dy 

%— 4 in n>0 WS (3.5). 
rt len fe 

oc (dp — ds) e o( Yen) 


Wir sahen bereits, daß für jede Stabstelle x die ursprüngliche Geschwindigkeit — v, durch die 
0) ar 
Stoßwelle nicht vollannulliert wird. Aus (3.5) folgt nun, daß En längs & = 0 und n = O positiv 
0 
ist. Da aber EL wegen (L,) auch eine Darstellung (L,) gestattet, ist dann überall el. 


Für jedes x sinkt also |v| monoton auf den Wert Null, und zwar für große t asymptotisch wie 
t-!12. Anders ist das Verhalten von o, das Bild 2 in dimensionsloser Darstellung durch Angabe 


c 
des o-Verlaufes für verschiedene t-Werte zeigt. Während nämlich für 5 x <2 der von der 


gestrichelt gezeichneten Stoßwelle aufgeprägte Druck mit wachsendem t monoton gegen — 0 


konvergiert, findet für m 2>2 zunächst eine Druckerhöhung statt, die physikalisch be- 
u 
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deutet, daß wegen der rascheren Spannungsrelaxation am Stabanfang schließlich eine Stauung 
der erzeugten Wellen auftritt. Für große t haben wir wieder asymptotisches Verhalten von 
— (0 + 0,) wie t-U2, 

Nachdem o bekannt ist, wird die bleibende Verformung e,; durch 


EL RL 
A u 


geliefert. Das Ergebnis der Rechnung zeigt Bild 3. 
Die gestrichelte Linie gibt dabei die Verformung für 


den Fall an, daß bei einem Stab der Länge |, 2 


die weitere Verformung durch die vom freien Ende 
zurücklaufende Entlastungswelle beendet wird. Die 
entsprechend verformte Gestalt eines zylindrischen 
Stabes nach der Reflexion sehen wir in Bild 4. 


Der von der festen Wand bis zur Zeit t über- 


t 
tragene Impuls I(t)= — [o(&=0)dt ergibt sich 
0 


aus (9.5) zu 


t 


PANETE 2 
I) =0-t+0e(n—n): |e ENAR (1552) 
g 2 


Für große t verhält sich das Integral wie 1/2. 


ur eınen Stab der Länge 
Lo ® 


feste wand 


REF Rare vor der Reflexion 


Bild 4. Bleibende Deformation eines zylindrischen Stabes endlicher Länge 


4. Spannung und Geschwindigkeit beim endlichen Stab 


Bei einem Stab endlicher Länge I, gelten die Ausführungen des vorigen Abschnittes so- 
lange, wie die vom freien Ende kommende Entlastungswelle AB von Bild 1 noch keinen Einfluß 
hat; d.h. (3.5) gilt in Gebiet II von Bild 1 einschließlich der Flanke & — &, — 0 der Entlastungs- 
welle AB, wobei 


C 
&=5, 2 Ye 1 dei 
gesetzt ist. Die Berechnung der Werte längs & = &, + 0 geschieht nun wieder durch Verwen- 
dung von (W) und (St) resp. längs und quer zu AB. Dabei wird.entsprechend der Spannungs- 


freiheit des freien Endes die Beziehung (W) längs &—= &,+ 0 zunächst als elastische angesetzt. 
Die einfache Rechnung führt auf 


o 3 
Fre TAT) | 
| für = 0 su. „u 


v 20, —v ; 
ee le 


HU U 
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Dabei ist zur Vereinfachung die Funktion 
. Kia EH y . N u: ” 
(8 y)= eV J,(2iYay) + J RA Nee ee 


eingeführt, von der man leicht zeigen kann?), daß sie bei festem zeine monoton mit y wachsende 
Funktion ist mit den speziellen Werten 


1 
9(,0)=e* und (a 9)=— .(1+e?rJ,(2i2) ..... (44). 
Nach (4.2) nimmt o längs &=&,+0 vom Werte Null beginnend monoton ab. Die Ent- 


lastung läßt also einen variablen Druck zurück, der zu neuen Wellen Anlaß gibt. Wenn o längs 
AB den Wert — o, unterschreitet, so wird (4. 2) Er Wir hätten dann eine auf AB be- 


ginnende, ins Gebiet I/II ver- 
EEE: ii 


laufende neue Grenze zwischen 


rein-elastischem und elasto- 
4 HH 
LEER SEEE 


hen, so daß geschlossene Aus- 
sagen nicht mehr möglich sind. 
Wir wollen uns daher auf den 
Fall beschränken, daß die Ent- 
lastungswelle bis zum Punkte B ‚ix 
von Bild 1in einen rein-elasti- 5 
schen Bereich führt, der natür- 

lich die inzwischen erreichten 
bleibenden Deformationen kon- 
serviert. Offen bleibt dabei 
vorläufig die Frage, ob trotz- 0 1 2 
dem vielleicht später der Stab 

wieder infolge der Wellenüber- ur 

lagerungen ganz oder teilweise 

in den plastischen Bereich gelangt. Die genannte Einschränkung ist wegen des monotonen 
Verhaltens von o längs AB gewährleistet, wenn |o(B)| <o, ist, was wegen (4.2) mit (4.4) 
führt auf: 


rakteristikenintegration gesche- | 
v0 


plastischem Verhalten. Die Ver- P 


folgung dieser Grenze kann nur 
ke 
3 


noch durch schrittweise Cha- 
en 


365) 


h —1 
Y<dgpi Mit a —=1+ & — ei + 5 e 2%), (238%) 12 (4.5). 
Der Verlauf von v,,.,/v, als Funktion der dimensionslosen Stablänge £, ist in Fig. 5 dargestellt. 
Verit/ vs besitzt ein Minimum vom Betrage 2,79 an der Stelle &, = 5,05 und wächst dann asym- 
ptotisch auf den Wert 3. Da bei v,>3  v, unser linearer Ansatz sowieso fraglich wird, ist auch 
insoweit unsere Beschränkung auf 1, < veris gerechtfertigt. 

Nachdem wir nun für das Gebiet I/II von Bild1 o und v längs AB kennen sowie die Rand- 
bedingung o (x = l,) = 0 haben, können wir o und vin III unter der stillschweigenden Vor- 
aussetzung des Verbleibens im elastischen Gebiete in elementarer Weise berechnen und er- 
halten: 

o 


Pc — %,) 
S EV 


%g gr UA %g at 


= 9 (Ed E — 80) — P (do: N) | 
HIOlllEr EN: ..1:(4.6). 


2) Es genügt nämlich y = e +" der Differentialgleichung Ya) = x und wird daher quellenmäßig darge- 
stellt durch ; 
u ei n 2 urF, 
eV = Il2ilEn) + SIo2iVE)ETAE+ SIEH ET E, 
0 0 


woraus sich bei& = n= zergibt: 
© a4 
2 (Il2ilat)e "FE=1-E°* I (2ir), 
0 


was zur Umformung von p(z, x) dient. 


Z. angew. Math. Mech. 


242 Richter, Elasto-plastische Reflexion eines Stabes Ba 98 Nr.7 Juli 1058 


Auf Grund der obengenannten Eigenschaften von g und wegen (4.5) zeigt (4.6), daß in III 
dauernd — 0,<o< O gilt, womit der rein-elastische Ansatz in diesem Gebiete gerechtfertigt 
ist. 

Um das anschließende Gebiet /V zu beherrschen, müssen wir wissen, ob nach Ankunft 
der Entlastungswelle in B der Stab sich von der Wand abhebt oder dort verbleibt. Da nämlich 
bei Bin III Druck herrscht, sagt der zugehörige Wert von v darüber noch nichts aus. Wir 
lassen daher ABin Breflektieren, um die neue Front BC zu erhalten; dies geschieht entweder 
unter der Bedingung v(2£= 0) —=0 bei Verbleiben oder der Bedingung o(&=0)=0 bei Ab- 
heben. In ZV muß dann im ersteren Falle o(B)< 0 und im zweiten Falle v(B)>0 sein. Die 
elementare Rechnung liefert diesbezüglich zunächst 

v<v mit = 1 All 2] rat 

8 

als Bedingung dafür, daß die Ankunft der Entlastungswelle am vorderen Stabende den Re- 
flexionsvorgang beendet. Es ist dies gleichzeitig die Bedingung dafür, daß sich das freie 
Ende nach der Entlastung von der Wand wegbewegt und gewissermaßen längs AB ein Abhebe- 
signal zum vorderen Ende sendet. Die Grenzgeschwindigkeit v; ist ebenfalls in Bild 5 mit ange- 
geben, wobei nur der allein interessierende Teil &, > 1,049 gezeichnet wurde, für den Y< verit 
ist. v, fällt monoton mit wachsendem £, bis zum Werte 2 v,. 


Es ergibt sich nun die folgende Fallunterscheidung: 


&%) %y< %ı. Diese Bedingung ist für &,< 1,049 stets erfüllt. Bei B findet Abheben des 
Stabes statt. Spannung und Geschwindigkeit in den Gebieten /V, V usw. ergeben sich im 
x-t-Diagramm wie üblich durch Spiegelung an x =1,/2 mit anschließender Translation in 
t-Richtung, wobei die Spannungen das Vorzeichen wechseln. Da Ill rein-elastisch war, bleibt 
der Stab im elastischen Zustand. Überdies überzeugt man sich leicht, daß gerade (4.7) garan- 
tiert, daß dauernd v(& = 0) > 0 ist, so daß der Stab auch später nicht nochmals infolge seiner 
Schwingungen die Wand berührt. 


ß) %<V<Verir. Dieser Fall kann nur für &,> 1,049 auftreten. Die neue Wellenfront 


BC, d-h. n = &,, ist hier unter der Voraussetzung v(x = 0) = 0 zu berechnen, wobei wir vor- 
behaltlich späterer Verifizierung die Rechnung zunächst rein elastisch durchführen. Man erhält 
dann in üblicher Weise für ound vin den Gebieten /V und V die Formeln: 


R 2%, —% 
RR En I ri Var 
» inIV <(4.8a) 
a: == P(Eo; 5) — (do, n— &0) 
0 8 
und 
o 
Eee ne AUT 2) be A | 
> RER? inV (4.8b). 
a en ae mar zn A ee 


Gemäß (4.8a) nimmt o von B bis D monoton zu mit o(B)<0 
und o(D)>0. Es gibt also eine bestimmte Stelle &=&, + 2, n=&, +2, 
wo an der Wand der Druck in Spannung übergeht und daher die 
Ablösung erfolgt. z bestimmt sich dabei aus der transzendenten 
Gleichung 


%— 2V% 
2 ‚er NT 4) 
9a, )= N —, (4.9) 
Bild 6 hängt also von &,und w/v, ab. Bei festem &, wächst zmonoton mit %.. 


An Stelle von Bild 1 verwenden wir nun das x-ct-Diagramm von Bild 6, in dem B’den 
Ablösepunkt kennzeichnet, von dem nunmehr eine Wirkungslinie ausgeht, oberhalb derer die 
gefundene Lösung für IV und V ungültig wird und durch eine andere zu ersetzen ist, bei der 
o(z = 0) = 0 als neue Randbedingung gilt. Über B’C’ hinweg sind o und v stetig fortzusetzen; 
längs B’C’ springen o und v natürlich am Schnittpunkt ? mit CD. Nehmen wir weiterhin zu- 
nächst an, daß das elastische Gebiet nicht mehr verlassen wird und daß die Randbedingung 
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o(©=0)=0 erhalten bleibt, so können wir elementar o und v zunächst in B’C’D’ und dann 
mit Hilfe des am Ende von «) genannten Spiegelungsprinzips in die weiteren anschließenden 
Dreiecke hinein berechnen, wobei gleichgelagerte Dreiecke gleiche Funktionswerte erhalten. 


Der Stab führt also nach der Ablösung periodische Schwingungen aus. Das Ergebnis dieser 
Rechnung ist: 


(4.8a) in BOPB'; (4.8b) in CC’ P, 


= _o(&, 6). len, uk) 


en) in B'PD, 
v 2, —v 
= pli E—E)+ Pl ne) + 
'o 8 U —d, 
= pl 2) pl + / 
ea 2 &0) — Pl&o» N — &0) sg | 2410), 
. in PC'OD, 
N = —9(&,E—28)+9(&,n— &,) 
0 8 
[03 
Er VERA EEE N zer) nn 
eco —v,) Pl, E — 2 5)+ Pd, n —2 80) 
2 ER EN DQOD 
5 ea) erlkon2)t 
di) 8 0 8 


Mit Hilfe der genannten Eigenschaften der Funktion @ und unter Beachtung von (4.5) 
überzeugt man sich leicht, daß in (4.10) dauernd |o|<< o, und v(z = 0) > 0 bleiben, so daß der 
Stab sich weder nochmals anlegt noch das elastische Gebiet verläßt. Damit sind diese beiden 
Voraussetzungen unserer Rechnung verifiziert. Es sei schließlich noch auf die bemerkenswerte 
Tatsache hingewiesen, daß das Abheben des Stabes von der Wand gerade in dem Zeitpunkt 
geschieht, wo an seinem vorderen Ende diejenige Wirkungslinie vom freien Ende ankommt, die 
dem Vorzeichenwechsel von v<0 zu v»> 0 entspricht. Wie im Falle «) gibt also auch hier das 
freie Ende das Signal zum Verlassen der Wand. 


In Bild 7 sind die Bewegungslinien des Stabes während und nach der Reflexion wiederge- 
geben, wobei Bild7anur die rein-elastische Reflexion als Vergleich bringt, während sich Bild 7 ba 
und 7bß auf die beiden geschilderten Möglichkeiten der elasto-plastischen Reflexion beziehen. 
Um die Unterschiede zeichnerisch deutlich zu machen, wurde dabei das Verhältnis v,/c = o,/E 
wesentlich überhöht. 


N 
il 


IR 


2; III .. NN 04 06 08 1012 1 161820 
0 02 04 06 08 1012 14 16 1820 ne 00 ed 12 19 16 182,0 welienfyonten 
5, x 7a u ame Wirkungslinien des Ablösevorganges 
77 An 
Bild 7a Bild 7b& Bild 7bß. Elasto-plastische Reflexion eines 


Stabes der Länge /,= 4u/oc bei 


Y<m< Vgrit 
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5. Reflexionsdauer und Restitutionskoeffizient 


Die Reflexionsdauer T ist im elastischen Falle gegeben durch T = 21,/c, oder in dimen- 
sionslosen Variablen 


00° e 
Au Z: & (3 ) 
Wir sahen bereits, daß (5.1) gültig bleibt bei v,< vy; insbesondere bei &,<{ 1,049. Ist dagegen 
Au 


U<m< Verie, so tritt für T eine Vergrößerung um —, - 2 ein, wo z durch (4.9) definiert ist. 
en 00 


Damit ergibt sich die in Bild 8 gezeigte Abhängigkeit der Reflexionsdauer von v, bei festge- 
haltener Stablänge £,- 


/ 


+ 
3“ = 
—{ 


> 
N 
& 
© 
NN 
Na 
AN 


Sr 
DERDRZEEE 
Esel] | 
ABERENETE 
VArAVARERTE 

ZEBANRRADAT 


ae 
SIEFIBENNE 
747 


eu 
Engarsssnano 


0 1 2 3 4 I Vol —— 
Bild 8. Reflexionsdauer Tin Abhängigkeit von o, Bild 9. Restitutionskoeffizient Rin ‘ 
bei fester Stablänge /, Abhängigkeit von v, bei fester Stablänge /, 


Der Restitutionskoeffizient R ist definiert als das Verhältnis R= v,/v, der absoluten 
Schwerpunktgeschwindigkeiten v, und v, vor und nach der Reflexion. Für ,< »,ist R=1. 
Bei v,> v, tritt jedoch ein Verlust an Translationsenergie auf, der die in Wärme verwandelte 
Verformungsarbeit und die gemäß Bild 7 nach der Reflexion sichtbare Schwingungsenergie 
deckt. Es muß also R< 1 werden. 


Nach dem Impulssatz ist nun v, bestimmt durch 
T 
uhe = — [o(a=0) dd — we; 
(0) 


wobei sich o(& = 0) für OS !< nr -&, aus (3.5) und —lalls 9, > v; ist — für - -&,<St<Taus 
(4.8a) ergibt. Die Rechnung liefert 


& vg 


+21) 2-42): (v0@. 


Dabei treten die Glieder der zweiten Zeile nur im Falle z>0, also bei v,>v; auf. Die Ab- 
hängigkeit des Restitutionskoeffizienten von v, bei fester Stablänge I, ist in Bild 9 dargestellt. 


Ri a (1-2): [neioea 
0 
0) 


Fingegangen am 9. Juli 1952. 
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KLEINE MITTEILUNGEN 


Ein einfaches graphisches Verfahren zur Be- 
stimmung des günstigsten Anstellwinkels für 
den Start von Flugzeugen. 


1. Ein startendes Flugzeug (vgl. Bild 1) rolle mit 
der Geschwindigkeit v» geradlinig auf einer ebenen 
Startbahn in x-Richtung. Der beschleunigende 
Motorenschub sei $, der Luftwiderstand W und der 
Rollwiderstand R. Senkrecht zur Bewegungsrich- 
tung wirken Auftrieb A und Fluggewicht @. Die 
Größen 8, W, R, A und G seien bekannt und mögen 


Bild 1 


in bekannter Weise vom Anstellwinkel & des Flug- 
zeugs abhängen. Dann genügt die Flugzeugbewegung 
nach dem Newtonschen Beschleunigungsgesetz 
der Gleichung (g = Erdbeschleunigung, i = Start- 


zeit) 
d(G 
—l— 9 = re 
& ») IS) W R 
oder 
G dv dG v 
ee — ne r 
FREE, Ss W-—-R an (+) 


Betrachtet man nun die rechte Seite von (+) als von 
der Beschleunigung dv/dt unabhängig !) (in Wirklich- 
keit wird z.B. der Rollwiderstand R einen Be- 
schleunigungsanteil enthalten), so kann man durch 
Variation des Anstellwinkels & des Flugzeugs ver- 
suchen, diesen Ausdruck, und damit auch die Be- 
schleunigung dv/dt, in jedem Augenblick möglichst 
groß zu machen, Damit erreicht man auch für jeden 
Zeitpunkt eine maximale Rollgeschwindigkeit, also 
auch kürzeste Startstrecke. 


Besonders einfach zu übersehen wird diese Aufgabe, 
wenn man für den Augenblick die Abhängigkeit des 
Schubs $ vom Anstellwinkel & vernachlässigt (die 
Verminderung des Schubs geht etwa mit cos a, hat 
also bei kleinen Anstellwinkeln & keine große Be- 


deutung). Da außerdem auch die Glieder und 


d 
—_ T = in (+) von & unabhängig sind, braucht man 


also nur — W— R zum Maximum, oder W-+-R 
zum Minimum zu machen. In der Flugmechanik 
wird nun meist an Stelle des Anstellwinkels & der 
Auftrieb A bzw. sein dimensionsloser Beiwert 


G= betrachtet (q=Staudruck, #—=Bezugsfläche). 
Dementsprechend sollen auch hier W-+-R, bzw. 


die Beiwerte cy-+ cr = IH über c„ statt über & 


gqF 
aufgetragen werden. Man erhält dann das Bild 2. 
Soll hierin c,„ + c, einen Minimalwert annehmen, so 
tritt dieser offenbar für dasjenige c„ = ca ein, in 
dem die nach rechts parallelverschobene (gestrichelte) 


1) Dies wird im folgenden die einzige Einschränkung bleiben, 
die getroffen werden muß. 


C,(Ca)-Kurve die Kurve cu(c.4) berührt. Infolge der 
eindeutigen Zuordnung von & zu diesen Kurven hat 
man damit auch den günstigsten Anstellwinkel 
x = 0a* gefunden. 

Aus dieser graphischen Darstellung läßt sich auch 
sofort ablesen, ob das betrachtete Problem eine, oder 
sogar mehrere Lösungen hat. In dem meist betrach- 
teten Sonderfall der Festreibung (Landflugzeuge), 
wenn also R = u(@ — A) ist (u = Reibungsbeiwert), 
wird die Kurve c, über c„ eine Gerade mit dem An- 


a — —-— 


£a 


Bild 2 


stieg 1/u und Bild 2 vereinfacht sich zu dem be- 
kannten ?) Bild 3. 

Es erschwert die Betrachtung nun nicht mehr 
wesentlich, wenn nachträglich auch noch die Kurve 


Bild 3 


des Schubs hinzugefügt wird, indem man etwa 
Ss 
qF 

2) Man vergleiche etwa Diemer, Verfahren zur graphischen 
Untersuchung der Startverhältnisse eines Flugzeugs. Zeitschrift 
für Flugtechnik und Motorluftschiffahrt, 17 (1926), 8.325 bis 
326. Doch ist bei Die mer — wie wohlin allen bisher bekannt- 
gewordenen Untersuchungen — die Richtigkeit dieser Kon- 
struktion nur für einen Spezialfall nachgewiesen, indem über 
die Abhängigkeit von 8, W, R, A und @ vom Anstellwinkel a, 


von der Geschwindigkeit» und der Startzeit t sehr einschrän- 
kende Annahmen getroffen werden. 


% + cz in Bild 2 statt c, allein aufträgt (« - 
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Da jedoch c, die in Bild 4 skizzierte Gestalt hat, wo 
die Lage c„ des Maximums allein von dem Winkel o 
zwischen Luftschraubenzugrichtung und Flugzeug- 
längsachse abhängt (siehe Bild 1), kann man auch o 
nachträglich so bestimmen, daßc, —= ca wird. Wegen 
E> —=0( für %, wird nämlich in diesem Falle die 

a 
Konstruktion von Bild 2 exakt richtig, auch bei Be- 
rücksichtigung des Schubs 8, und außerdem hat man 
dann den kleinstmöglichen Startweg gefunden, falls 
die beiden Parameter x und o frei verfügbar 
sind 3). 

Rückblickend sei ausdrücklich festgestellt, daß 
für das geschilderte Verfahren keine besonderen An- 
nahmen bezüglich der verschiedenen auftretenden 


Bild 4 


Kräfte notwendig sind. Es gilt sowohl für Wasser- 
als auch für Landflugzeuge, für Luftschrauben- 
antrieb oder Schuberzeugung durch Strablapparate, 
unter Berücksichtigung des Gewichtsverlusts durch 
Brennstoffverbrauch und bei beliebiger Abhängigkeit 
der Luft-, Reibungs- und Schubkräfte von Startweg, 
Rollgeschwindigkeit und Zeit; wenn nur die rechte 
Seite der Gleichung (+) als von der Beschleunigung 
unabhängig angenommen wird. Das Verfahren 
bleibt also immer gültig, ob über eine bestimmte Zeit 
ein zusätzlicher Startschub durch Hilfsraketen er- 
zeugt wird, ob bei einer gewissen Rollgeschwindigkeit 
eine Startklappe ausgefahren wird ®), oder ob nach 
einer bekannten Rollstrecke eine Betonbahn zu Ende 
ist und das Flugzeug auf einer Grasnarbe mit größe- 
rer Reibung weiterrollt. Ebensowenig ändern sich 
diese Überlegungen durch den Einfluß von Boden- 
nähe und Reynoldsscher Zahl auf Widerstand 
und Auftrieb, durch Wind oder eine Neigung der 
Startbahn, ja der Startweg braucht noch nicht ein- 
mal geradlinig zu verlaufen, sondern kann beliebige 
Kurven beschreiben. 

2. Nachdem so die optimale Flugzeuganstellung 
während des ersten Teils des Starts gefunden ist, er- 
hebt sich die Frage nach den günstigsten Abhebe- 
bedingungen. Ist es z. B. zweckmäßig, so früh wie 
möglich abzuheben 5), oder soll man etwa bis zur 


®») Da die Wahldes Winkels o im allgemeinen auch noch durch 
andere Gesichtspunkte eingeengt wird, braucht es allerdings oft 
nicht möglich zu sein, gerade den günstigsten Winkel o* einzu- 
bauen. Wegen der geringen Veränderlichkeit des Schubs im 
interessierenden o-Bereich ist die dadurch hervorgerufene Start- 
verlängerung allerdings meist vernachlässigbar klein. 

*) Nicht behandelt wird hier die Frage, wann ein solcher 
Zusatzschub am günstigsten einzusetzen habe, oder bei welcher 
Geschwindigkeit man am zweckmäßigsten Starthilfsklappen 
ausfährt. 

®) Diese Forderung wird z. B. (ohne Begründung) aufgestellt 
bei A. Pröll, Die Startstrecke bei Flugzeugen. Zeitschrift 
für Flugtechnik und Motorluftschiffahrt, 17 (1926), S.316 
bis 322. 


Geschwindigkeit des besten Steigens rollen, um erst 
dann in den Steigflug überzugehen °)? Schon eine 
kurze Überlegung zeigt, daß diese allgemeine Frage- 
stellung sicb nicht ohne weitere Verabredungen 
lösen läßt. Nimmt man nämlich z. B. ein räumlich 
begrenztes Startfeld an, so ist es offenbar notwendig, 
spätestens am Flugplatzende abzuheben, wenn der 
Start überhaupt möglich sein soll. Durch solche 
äußeren Umstände kann also erzwungen werden, daß 
der zweckmäßigste Abhebepunkt an einer anderen 
Stelle liegt, als sich aus irgendeiner Optimums- 
betrachtung ergibt. 

Um trotz dieser grundsätzlichen Schwierigkeiten 
eine Problemstellung zu finden, die einigermaßen 
mit der Wirklichkeit übereinstimmt, soll angenommen 


Bild 5 


werden, daß das Flugzeug bis zu dem Punkte A 
(vgl. Bild 5) rollt, dort abhebt und zunächst im 
Horizontalflug dicht über der Erdoberfläche seine 
Geschwindigkeit steigert, um erst ab Punkt Bin den 
Steigflug überzugehen. Eine einfache Betrachtung 
an Hand von Bild 2 wird dann zeigen, daß sich der 
Punkt A automatisch allein aus der Forderung be- 
stimmt, während des Rollens immer möglichst 
günstig angestellt zu sein, oder, anders ausgedrückt, 
daß im Punkt A günstigster Roll- und Abhebe- 
anstellwinkel übereinstimmen. Für den Flugzeug- 
führer bedeutet das die einfache Forderung, das 
Flugzeug ‚von selbst“ abbeben zu lassen und es 
nicht etwa frühzeitig vom Boden wegzureißen ?). 


Bild 6 


Man erkennt diese Behauptung an Bild 2 und 6, 
wenn man bedenkt, daß längs der Strecke A — B 
von Bild5 die Gleichung (+) unverändert die Be- 
wegung des Flugzeugs beschreibt. Die Polare c„(c«) 
und die Kurve c,(cz) des Flugzeugs können sich näm- 
lich bei wachsender Geschwindigkeit v höchstens 
stetig verändern. Dabei rückt der Schnittpunkt P 
der Reibungskurve c,(c,) mit der Ordinatenachse 


°%) Man vergleiche Diemer [siehe Fußnote ?)], der dieses 
Vorgehen als zweckmäßig empfiehlt. 

?) Das steht nicht in Widerspruch zu der Erfahrungstatsache, 
daß schwerbeladene Seeflugzeuge u. U. erst durch periodische 
Anstellwinkelveränderungen ‚‚aufgeschaukelt‘‘ werden müssen, 
um überhaupt abheben zu können. Durch solch eine Maßnahme 
kann in divsem Falle die Reibungskraft R in ihrer Größe und in 
ihrem Charakter günstig verändert werden; was hier jedoch 
nieht zur Debatte steht. 
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stetig nach unten, um schließlich beim günstigsten 
Abhebe-c, auf den Berührungspunkt der verschobe- 
nen (gestrichelten) c,(c,)-Kurve mit der Polare cy(ca) 
zu treffen. Wie man aus dieser Darstellung ent- 
nimmt, ist also bei dem nach Bild 2 gewonnenen c* 
in jedem Augenblick des Starts — einschließlich des 
Abhebevorgangs — die Beschleunigung am größten, 
also auch die benötigte Startstrecke am kleinsten. 
3. Es ist nun sehr wohl denkbar, daß das Flugzeug 
durch äußere Umstände gezwungen ist, während des 
Anrollvorgangs nicht mit dem günstigsten Anstell- 
winkel nach Bild 2 zu starten. Durch einen Blick 
auf Bild 21äßt sich sofort wieder erkennen, wie dann 
der günstigste Abhebevorgang verläuft (Bild 7): Sei 


Bild 7 


z.B. während des Rollens stets 4 = iu, < ca Vvor- 
geschrieben, so wird man im Abhebepunkt A den 
Anstellwinkel & möglichst sprunghaft erhöhen, bis 
Cq = Ca, erreicht ist und das Flugzeug fliegt. Dann 
ist offenbar kein Verlust an Beschleunigung, und 
also auch keine Startstreckenverlängerung einge- 
treten, und bei dem anschließenden Fahrtaufholen 
im Waagerechtflug steigt nach Bild 7 die Beschleuni- 
gung i. a. stetig noch weiter an. 


Ist dagegen während des Rollens etwa 
* 
Ca = Ca, > Ca 


vorgeschrieben, so wird man wieder, wie bei 2., den 
Anstellwinkel zum Abheben nicht erhöhen. 


Es ergeben sich also je nachdem zwei ganz ver- 
schiedene Richtlinien für das Verhalten des Flugzeug- 
führers: Liegt das Roll-c, über dem günstigsten 
Wert ca, oder erreicht es ihn höchstens, so ist es 
zweckmäßig, das Flugzeug auch mit diesem Wert 
abheben zu lassen ; liegt es dagegen darunter, so wird 
zum Start der Anstellwinkel am besten ruckartig 
erhöht. Zu welchem Zeitpunkt dies zu geschehen hat 
und um wieviel gezogen werden muß, geht aus Bild 7 
hervor. 

4. Beispiel: Nähert man die Polare, wie dies oft 
geschieht, durch eine Parabel 


k 
= + =;ä cH an 
(A = b?/F = Flügelstreckung mit b = Flügelspann- 
weite und # = Flügelfläche, der Faktor k gibt die Ab- 
weichung gegenüber der elliptischen Auftriebsvertei- 
lung an), und ersetzt man c, + c, durch eine Gerade 


+ta=Ü+ua, 


so wird 
nn 


“72% 


aus Bild 7 liegen dann symmetrisch zu c*, 


us. Ca, und 


Ca, 


d.h. es ist 
k . . ı%* 
Cz — Ca, = Cay _— cz P 


Für elliptische Flügel gilt als Zusammenhang 
zwischen c„ und & näherungsweise, 


= 2un, rn 5 n ü, 


wo n der Profilwirkungsgrad ist. Damit erhält man 
schließlich als günstigsten Roll- und Abhebeanstell- 
winkel unter den geschilderten speziellen Annahmen 


A-+2n 
% 
5 4kn at; 


Zahlenbeispiel: Für k=en=1, A=1W, 
findet man 


u = 1/30 


o* = 0,12 5,73°, 


Cördoba (Argentinien). Karl Nickel. 


Elementarer Beweis für die Partialbruch- 
zerlegung des Cotangens. 

Wir geben im folgenden eine vereinfachte Beweis- 
anordnung für diese Zerlegung: der erste Teil ist der 
übliche, und stimmt z. B. mit der Darstellung bei 
Knopp!) überein; im zweiten Teil wird die Über- 
einstimmung der Reihe mit dem Cotangens jedoch 
in anderer Weise erschlossen. 

1. Es ist 
„608 20 


5 cos?x — sin? 
"sin?r "2sinz-cos® 


= ctgae —tgr 


(1) 


2ctg2r7 = 


und wegen 


also 


Mithin genügt die Funktion ctg ar x mit der Periode 1 
der Funktionalgleichung 


ES DZ RICH 2) 
cetgn © wir cy FH etg 3 ir ; 


Wendet man diese Funktionalgleichbung mit dem 
Pluszeichen auf jedes Glied rechts an, so folgt 


. (3). 


1 2 (2-1 
cetgra = ach --ctg Di, a ) | 
4) 
n(2 -- 2) (a -+ 3) ( 
N PX +ctg — 52 | 
und allgemein 
71 


l 
ctg na = gan 
er ; 1 
oder, wenn man für die Glieder von » = 2" — lab 
einschließlich das Argument z durch z — 1 ersetzt, 
d.h. » durch » — 2%, und symmetrisch zusammen- 
faßt 


gn—1_] 


1 nw a(® +») 
ganz Ag an a: 1; 10 
Br (6). 
n(2 —v) IH 


1) K. Knopp, Theorie und Anwendung der unendlichen 
Reihen, 4. Auflage, Berlin 1947, 8. 209—212. 
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Führt man rechts den Grenzübergang n — © glied- 
weise durch, unbeschadet der Frage, ob dies erlaubt 
sei, so erhält man 


Int ST 18, - 
ee, ER T), 


wor #0, +l, +2,...sei. Diese Reihe konvergiert 
für jedes solche x, und sogar gleichmäßig in jedem 
beschränkten abgeschlossenen Gebiet, das jene Punkte 
nicht enthält; ihre Summe sei C(z). Dann ist C(«) 
an jeder Stelle 40, +1, +2,...stetig. Zu zeigen 
ist 


cdgne=(Ü(e). (8). 
2. Dieser Nachweis erfolgt nun so: 


Man zeigt sofort, daß C(x) genau wie ctgn x die 
Periode 1 hat und derselben Frnktionalgleichung (3) 


genügt: 
ae 


Mithin gilt dasselbe für die Differenz 
h(a) =ctgne —((e): 
h(@-+1)=hl(e) 


ee) 


Rz) Qu 2, > 0er 12), 
und also infolge der Periodizität von h (x) auch 
h(z2)>0 für z>1—0 Ela). 
In0o<e<e<s1-— eist h(x) stetig, und also auf 
Grund von (12) und (13) zu einer n0O<z<<1 ste- 
tigen Funktion h(x) ergänzbar. Infolge der Periodizi- 
tät ist damit h(x) als stetige Funktion für alle x ge- 
geben; speziell ist A(x) integrierbar. Dieselbe Über- 
legung, die von (3) zu (5) führte, ergibt, angewandt 
auf irgendein x, aus (11) 
1 
z-+v 
(3) 


 < 
h (®) —yon 
Für n— w ergibt sich gemäß der Definition des 
Integrals 


ferner gilt 


. (14). 


1 
h\) la) dt. d.h. h(2) = consiralb) 


und, da A(0) = 0, mithin h(2) = 0, ctgne = C(k) 
oder 


I gt 1 
TE SUE TEN a az a (16). 


Ernst Mohr. 


Berlin-Charlottenburg. 


Zur graphischen Berechnung von Polynomen. 

Das Verfahren vonH.Behmann!) kommt mit 
einem Minimum an Hilfslinien aus, die gar nicht ge- 
zeichnet zu werden brauchen, sondern deren Schnitt- 
punkte bloß mit einer dem gewählten Argument ent- 
sprechenden Gerade zu markieren sind. Man kann 
aber eine schon von A. Winckler?) aus dem 
Mesolabium des Eratosthenes hergeleitete 
Konstruktion, die vor einigen Jahren F. W. Palm?) 
eingehend untersucht hat, durch Verwendung von 
Pauspapier noch so vereinfachen, daß man auf 
diesem für einen gegebenen Argumentwert nureine 
einzige Gerade zu ziehen hat, sobald ein für 


alle Fälle brauchbares Strahlbüschel und eine 
Parallenschar gezeichnet vorliegen, die auf eine mit 
einer regulären Skala versehene Achse senkrecht 
steht. 

Um 


fe)=%' 2" a, zn l + +4 in © + An 

zu berechnen, bedient man sich des bekannten Sche- 
mas: 
Ay a, a Fol An, 

Go © Pı TV’ Pn- %' Pn-ı 

2122.05 ... Pnıı MmM=fl®), 
wobei die p7; auf die aus der Abbildung ersichtliche 
Weise graphisch gebildet werden. Es ist dort 
OE,=0E,=1, ferner mit Berücksichtigung des 
Vorzeichens 4,0, = ar, also 00, =p, und für 
x = E„P soll f(x) (n = 4) gefunden werden: f(2)=00,. 


F@) = 16x*+ 0,8x°-14x2-0,4x+0,6 
f (0,4) =0,31 


Zur praktischen Durchführung zieht man auf Pavs- 
papier eine Gerade »/, trägt auf ihr von einem Punkt 
O0’ aus mit der Längeneinheit 0OE, die Koeffizienten 
a; ab, und zwar die positiven nach unten und die 
negativennach oben. Die Endpunkte dieser Strecken 
seien A,, A], 43.... Das Pauspapier legt man dann 
so auf das Blatt mit dem Büschel, dessen Strahlen 
nach x beziffert sind, daß die Gerade n’ mit der 
n-Achse und der Punkt 0’ mit E, zusammenfallen, 
und ziebt durch 0’ die Gerade «’ nach dem Schnitt- 
punkt des Strahles x mit der n-Parallelen durch E,. 
Nun verschiebt man das Pauspapier so, daß n’ stets 
über n bleibt und zunächst A/ sich mit O deckt. 
Mittels der &-Parallenschar sucht man -— eventuell 
durch Interpolation nach dem Augenmaß — den auf 
gleicher Höhe wie P, = (2, x’) liegenden Punkt A, 
der n-Skala, den man sich aber nur für den folgenden 
Schritt zu merken braucht. Auf A, legt man dann 
A7, merkt sich wieder den auf gleicher Höhe wie 
P, = (x, x’) liegenden Punkt A, der n-Skala, auf den 
man 4% legt usw. Die Ordinate des Punktes 4}, 
liefert schließlich beim letzten Schritt den gesuchten 
Wert von f(x). 


1) Zur graphischen Behandlung der ganzen rationalen Funk- 
tion. Z. angew. Math. Mech. 11 (1931), 8. 463. 

2) Geometrische Konstruktion rationaler Polynome. Sitz: 
Ber. d. Wiener Akad. d. Wiss. Ila, 53 (1866), S. 326. 

3) Anwendung und Verallgemeinerung des graphischen Ver- 
fahrens von Winckler. Sitz.- Ber. d. Österr. Akad. d. Wiss. IIa 
157, (1949), S. 275. 
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Wegen 
1 = 
en (4) =Ay20 HA, an 1-4... +a,2-+0, = f(x) 


braucht man das universale Strahlbüschel nur für 
—1<sSx<s1zu zeichnen und hat beim Übergang 
von x zu 1/2 bloß die Reihenfolge der Punkte A, auf 
dem Pauspapier umzukehren. 

Für variables x beschreibt der Punkt P(£,n) eine 
Kurve mit der Parameterdarstellung 


1 


g=—f(a); 


n=f(?), 


aus der folgt: = f Bi) Sie ist also rational von der 
Ss 


(r + 1)-ten Ordnung und hat im Ursprung einen 
n-fachen Punkt mit den Tangenten 7 = 7, : £&, wobei 
He .0: 


Wien. 


A. Huber. 


Zur Düsenströmung mit Reibung. 


Bei der thermodynamischen Berechnung von 
Strömungsmaschinen wird i.a. die Reibung an den 
Kanalwänden durch Einführung eines sog. Ge- 
schwindigkeitsbeiwertes unter der Annahme berück- 
sichtigt, daß bei konstantem Druck im Spalt zwischen 
Düse und Laufrad eine dem Geschwindigkeitsverlust 
entsprechende Enthalpiezunahme und Entropie- 
vermehrung des Arbeitsmittels eintritt; der Be- 
stimmung des Durchsatzes werden üblicherweise 
unter Verwendung eines sog. Sicherheitsfaktors die- 
jenigen Zustandsgrößen im engsten Querschnitt 
zugrunde gelegt, die sich bei isentropischer Strömung 
unmittelbar aus den kritischen Werten ergeben. — 
In der vorliegenden kurzen Notiz ist versucht worden, 
nach Ableitung der für die Düsenströmung mit Rei- 
bung gültigen Beziehungen eine für technische 
Zwecke geeignete Näherungslösung aufzustellen. 

Bei der hier vorausgesetzten Expansionsströmung 
möge die an der Kanalwand durch Reibung erzeugte 
Wärme dem Gas verlustlos wieder zugeführt werden; 
der Kanal möge die Gestalt einer Lavaldüse haben, 
das Druckgefälle zwischen Düseneintritt und -aus- 
tritt sei als überkritisch vorausgesetzt. Der gas- 
dynamische Zustand im engsten Querschnitt und 
somit auch der Massendurchsatz ist infolge der 
adiabatischen pichtisentropischen Strömung zu- 
nächst unbekannt; die Bestimmung dieser Werte 
kann nicht durch die Lösung einer einfachen Extre- 
malaufgabe erfolgen, wie dies bei der Betrachtung 
der reibungsfreien Strömung geschieht. 

Wird für die von der Düsenoberfläche dO = 
nD-dx (kreisförmiger Querschnitt, D: Düsen- 
durchmesser, x: Düsenabzisse, in Strömungsrichtung 
positiv gerechnet) herrührende Reibungskraft dR der 
Ansatz 


dR = er DS w- de 
(er: Reibungsbeiwert, o: Dichte, w: Geschwindigkeit) 
gemacht, so gilt wegen dm = oF - dx (m: Masse des 


Mediums, F: Düsenquerschnitt) für die der Massen- 
einheit durch die Reibung zugeführte Energie de: 


dR TC D 
mE en .de: 
e=7 da RZWT dx; 


daher ist nach dem 2. Hauptsatz der Wärmelehre 
Be 
T ds ROT de. ERST) 


(T: Temperatur, s: Entropie je Masseneinheit). In 
der Kontinvitätsbeziehung 


Fow= m(Tmin, Smin) =&' m(cr = 0) - (2) 


soll die in der Zeiteinheit durch einen Kanalquer- 
schnitt strömende Masse ri, die bei gegebenem Aus- 
gangszustand vom Zustand im engsten Düsen- 
querschnitt (Index ‚‚min“) abhängt, durch den Ver- 
lustfaktor £ auf die leicht berechenbare isentropische 
Strömung bezogen werden. Die thermische Zustands- 
gleichung, in den Variablen og, T und s geschrieben 


1 
uns A en N ile0y 
Fun T, E ER 


(Cp, v: spezifische Wärmen je Masseneinheit, der In- 
dex „0“ bezieht sich auf den Ausgangszustand) und 
die Energiebeziehung 


w“— u =2%'(7,—T) 


. (3) 


= (4), 


vervollständigen das zur Bestimmung der Unbe- 

kannten 7, o, s und w notwendige Gleichungssystem. 
Die Düsenabmessungen sollen auf den engsten 

Querschnitt bezogen werden. - 


Ö(E) == D(E)/.Dmin . 


Ferner soll die sog. Kesseltemperatur (Ruhezustand) 
eingeführt 


E = 2].Dmin , 


Too = To + wil2 cy 
und 


® = T/Too; = (8 — So)/(Cp — ©) 


gesetzt werden. Dann nehmen die aus (1) und (4) 

bzw. (2) und (4) sich ergebenden Beziehungen 
die Form an: 

4x 

”»—|1 


1-9 1 do 


TE 
1 
Kr =-d’ 11 —-92.8,...X6) 


mit 


1 Bi 1 

a 2 \x— 11/2 — 1[00/ To\* 1 n 
&= (4) Ver = 7) ) 9 
Nach (5) und (6) gehorcht die Temperatur im Kanal 


der für einen vorgegebenen Düsenquerschnitt nu- 
merisch zu lösenden Gleichung 


Ar Es. dd 
2 1-9 x —1 d)deE 


6 4x 1 


- (8) 


in welcher der Abwertungsfaktor & nicht auftritt 
(6° = dö/d£). Aus (8) ergeben sich nachstehende 
Folgerungen: 

a) Da die örtliche Schallgeschwindigkeit c lediglich 
von der Temperatur und die Strömungsgeschwindig- 
keit nach (4) nur vom Temperaturgefälle abhängt, 
wird die kritische Temperatur 


Or = - (9) 


“+1 
durch die Reibung nicht beeinflußt; daher wird 
nach (8) der kritische Zustand (c = w) hinter dem 


engsten Düsenquerschnitt — in Strömungsrichtung 
gerechnet — an der Stelle &7, erreicht, an der 

Ö=xcr (10) 
bzw. 


d2ö 


1 
Err 7 Emin + 27 KER (=) ' | 
mın 


; . (11) 
= Emin + rk "min 
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ist, wobei "min = 2 Rmin: Dmin gesetzt ist und Amin 
den Krümmungshalbmesser am engsten Düsen- 
querschnitt bedeutet. 

b) Für die Temperatur din im engsten Düsen- 
querschnitt (6° = 0, ö = 1) liefert der Ansatz 


H do 
Dr = din + (kr — a) 
min 
unter Berücksichtigung von (8), (9), (10) und (11) 


nach umfangreichen elementaren Umformungen den 
Ausdruck 


1 + ren rn | 


Omin = 


. (12). 


re 

1+xecr V "min 
Nach Ermittlung der Temperavurverteilung nach (8) 
kann die Entropie nach (5) durch einfache Quadratur 
bestimmt werden, wobei sich der thermodynamische 
Zustand des Mediums im engsten Querschnitt und 
nach (6) der Verlustfaktor £, also der Gasdurchsatz 
ergibt. 

Zur Vermeidung der ausgedehnten numerischen 
Rschnungen kann unter ersichtlichen Vernach- 
lässigungen ((o = 0 in (6)) ein hinreichend befrie- 
digender Näherungswert der Entropie berechnet 
werden, zu welchem Zweck ein Ersatz der in (6) auf- 
tretenden Temperaturfunktion notwendig ist: 

Im Anschluß an Stodola kann im unterkriti- 
schen bzw. schwachüberkritischen Gebiet die für den 
Durchsatz charakteristische Funktion 


Der 


1 
fy=9n!/1-9 


durch eine Ellipsengleichung mit guter Annäherung 
wiedergegeben werden: 


20-1 - Be) 


wobei die Konstante fz, mit der in (7) eingeführten 
Konstanten K wie folgt zusammenhängt: 


(1a) 


1 


1 
=1 dh —ı 
ee Tee 
hr kr Y kr GR 
Mit (13) liefert (6) und (5) einen ersten Näherungswert 
für die Entropie im engsten Querschnitt: 


Smin 


Omin=4xcrR 


mit 


Für technische Berechnungen dürfte bereits der aus 
(14) in Verbindung mit (12) bestimmte thermo- 
dynamische Zustand im engsten Kanalquerschnitt 
und der damit nach (6) und (2) sich ergebende Gas- 
durchsatz hinreichend genau ermittelt sein. Eine 
weitere Verbesserung könnte bei variabler oberer 
Grenze des Integrals in (14) wie üblich nach der 
Methode der sukzessiven Approximation gewonnen 
werden; die hierzu notwendige Rechenarbeit würde 
jedoch an Umfang nicht derjenigen nachstehen, die 
bei der ohne Vernachlässigung durchgeführten Rech- 
nung erforderlich ist. 


Bei einer nicht erweiterten Düse bzw. einer Laval- 
düse mit unterkritischem Gegendruck tritt an Stelle 
von (12) die Bedingung, daß am Kanalaustritt der 
vorgegebene Spaltdruck herrschen muß; in Ver- 
bindung mit (14) ergibt sich aus dieser Forderung 
eindeutig der gasdynamische Zustand am Düsen- 
austritt. 


Berlin-Lichtenberg. W. Bader. 


Differentiation durch Integration 

Zur Differentiation empirischer als Tabelle gegebener 
Funktionen werden in der Regel die Ableitungen der 
Interpolationspolynome benutzt. Meßfehler der em- 
pirischen Funktion lassen dabei die Differential- 
quotienten nur ungenau angeben. Verbindet man je- 
doch die Differentiation mit einer Integration, so werden 
durch diese die Meßfehler ausgeglichen. 

Für die Ableitung y’ der gegebenen Funktion % setzt 
man ein Polynom b, +52 +: + b„xz?% an und be- 
stimmt dessen Beiwerte b aus dem Minimalprinzip 


+1 
V’=/fSt{b,+bs+*--+bnar — y’}?®de = Min. 
1 


Durch Differentiation von V nach den b, erhält man 
dien + l linearen Gleichungen für die bo: 


a 
Sa? bu +bhr +: + bnen)de = 
—1 
+1 +1 
= [ze y]) —o-fatTt-yde=Ko 
— 1 
005 12): 


Für n = 3 ergeben sich z. B. die Beiwerte 


I 


3 15 7 
by 3 ®@R-5K); = BKL =TR) 


15 35 
=, 0 — 0); b=3(5K—-3Kı). 


Nach Berechnung der Beiwerte bo ermöglicht der Ver- 
P7 

gleich von / (bu, +b, 2 +: + b„ar)-de mit y(x) 

eine Fehlerabschätzung. 


Hannover. H.v. Sanden. 


BUCHBESPRECHUNGEN 


M. Bergsträsser, Evolventengeometrie 
für Stirnrädergetriebe. (VDI-Forschungs- 
heft 436.) 24 + 38 Seiten m. 41 Abb. u. 3 Zablen- 
tafeln. Düsseldorf 1952. Deutscher Ingenieur-Verlag 
G.m.b.H. Preis 15,— DM, für VDI-Mitglieder 
13,50 DM. 

Das Heft gibt eine klare Gliederung der Grund- 
begriffe und Paarungsgleichungen für Geradzahn- 
Stirnräder, der Vergleichsbeiwerte, des Berechnungs- 
ganges und der Verzahnungsbedingung, der allgemei- 


nen Korrektur schrägverzabnter Stirnräder und der 
mathematischen Grundlagen. Die Untersuchung über 
die Geometrie der Verzahnungen wird der Fachwelt 
aus dem Grunde vorgelegt, weil in ihr eine von der 
üblichen Darstellungsweise abweichende Frage- 
stellung erörtert wird, die zu nützlichen Paarungs- 
gleichungen für die Behandlung der allgemeinen 
Evolventenverzahnung führt. Über den Wert der Kor- 
rektur der Verzahnung ist man sich in der Praxis seit 
langem im klaren. Es ist zu begrüßen, daß in dem vor- 
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liegenden Heft alle geometrischen und rechnerischen 
Grundlagen in vorbildlicher Weise zusammengestellt 
sind und damit dem Fachmann zugänglich gemacht 
werden. Der günstige Einfluß, der durch die Profil- 
verschiebung bei Gradzahnrädern hinsichtlich einer 
Verbesserung der Biegefestigkeit der Zähne, einer Ver- 
ringerung der Walzenpressung, einer Verbesserung 
der Laufeigenschaften sowie der Einhaltung eines ge- 
gebenen Achsabstandes bei Verwendung genormter 
Werkzeuge erreicht werden kann, wird mit hinreichen- 
der Deutlichkeit an mehreren Beispielen hervor- 
gehoben, welche gleichzeitig zeigen, wie das allgemeine 
Korrekturverfahren auch bei Schrägzahn-Stirnrädern 
zu besseren Zabnformen führt als die Null-Verzahnung. 

Hervorzuheben ist, daß der Verfasser — soweit 
möglich — auf die DIN-Normen der Zahnräder Bezug 
nimmt und in den verschiedenen Abschnitten durch 
kritische Betrachtung der verschiedenen Möglich- 
keiten den Leser bei der sinnvollen Verarbeitung der 
durch die Zahnformenkorrektur gegebenen Vorteile 
unterstützt. 


Dresden. Lichtenheldt. 


Prof. Dr. A. Adler, Wahrscheinlichkeit 
und Wissenschaft. 29S. m. 2Abb. Bern 
1953. Verlag Paul Haupt. Preis kart. 2,10 Schw. Fr. 

Das Heftchen enthält eine Berner Rektoratsrede. 
In ihr wird von den verschiedenen Definitionen der 
Wahrscheinlichkeit zunächst die klassische von 
Laplace, dann die Häufigkeitsdefinition von 
v.Mises unddievon Cram&r behandelt, weiter 
einige Verteilungen wie die Binomial-, die Normal- 
und die Poissonsche Verteilung und ihre Bedeutung 
besprochen und schließlich der x?-Test und seine An- 
wendung unter anderem auf die schweizerischen 
Sterbetafeln erläutert. In dem klar geschriebenen gut 
lesbaren Heft wird vor allem die immer mehr wach- 
sende Bedeutung der Wahrscheinlichkeitsrechnung 
für andere Wissenschaftsgebiete betont. 


Dresden. Willers. 


Jenaer Jahrbuch 1952 (Wissenschaftliche Ver- 
öffentlichungen des Zeißwerkes). XI + 
266 S. m. 173 Abb. Jena 1952. Kommissionsverlag 
Gustav Fischer. Preis geb. 18,— DM 

Zum dritten Mal erscheint das Jenaer Jabrbuch, 
dessen wissenschaftlicher Teil dieses Mal durch eine 
Darstellung der geschichtlichen Entwicklung mikro- 
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pbotographischer Apparate im Zeißwerk eingeleitet 
wird. Außerdem entbält es u. a. Arbeiten aus den Ge- 
bieten dünner Schichten, der angewandten Optik, der 
lichtelektrischen Erscheinungen und der Zehntelmilli- 
meterwellen. Den Leser dieser Zeitschrift werden ins- 
besondere die Arbeiten aus dem Gebiet der angewandten 
Mathematik und Mechanik interessieren, so die Arbeiten 
vonKneschke: Über den Zusammenhang differen- 
tieller Integrationsprobleme mit den Interpolations- 
theoremen, von Hasselmeier: Die Kurzkegel- 
zentrierung und ihre Bemessung und vonBauer und 
Nessler: Wärmedurchgang durch anisotropePlatten. 


Dresden. Willers. 


E. v. Angerer, Technische Kunstgriffe 
beiphysikalischen Untersuchungen. 
8., vollständig neubearbeitete und erweiterte Auflage. 
Herausgegeben von H. Ebert, Braunschweig, 
unter Mitarbeit zahlreicher Fachwissenschaftler. 
352 S. mit 118 Abb. Braunschweig 1952, Verlag Fr. 
Vieweg & Sohn. Preis geb. 14,80 DM. 


Die neue Auflage des wohlbekannten Buches von 
Angerer ist vollständig umgearbeitet worden. Aus 
dem Werke, das aus der Feder eines einzigen Autors 
stammt, ist jetzt eine Gemeinschaftsarbeit geworden, 
an der fast 30 Mitarbeiter tätig gewesen sind. Diese 
haben jeweils kleinere Abschnitte verfaßt. Bei dem 
Anspruch des Buches, Kunstgriffe aus dem ganzen 
Gebiet der Physik zu bringen, drängte sich diese Ar- 
beitsteilung natürlich dem Herausgeber auf. Auf- 
gabenstellung und Charakter des Werkes sind trotz- 
dem so geblieben, wie sie durch Angerer vorge- 
zeichnet waren. Sorgfältige Darstellung der einzelnen 
Kunstgriffe zeichnen das Buch aus; jeder Kunstgriff 
sollte möglichst von dem Autor selbst erprobt sein. 
Umfangreiche Literaturangaben weisen auf ausführ- 
liche Arbeiten bin. Ein ausführliches Sachverzeichnis 
und ein Bezugsquellennachweis sind vorhanden. 
Durch verhältnismäßig kleinen Druck und durch 
Dünnpapier hat man erreicht, daß trotz des umfang- 
reichen Stoffes die äußere Gestalt eines Taschenbuches 
bewahrt geblieben ist. Jedem Physiker ist geläufig, 
wie gut das Buch von Angerer hinsichtlich Um- 
fang des gebotenen Stoffes und Zuverlässigkeit der 
Angaben ist. Diese Vorzüge hat es auch in der neuen 
Auflage behalten. 


Dresden. A. Recknagel. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


# 
Bei der Schriftleitung sind folgende Bücher eingegangen, ausführliche Besprechung bleibt vorbehalten. 


H. Rutishauser, Automatische Rechen- 
plananfertigung bei programmge- 
steuerten Rechenmaschinen (Mitteilun- 
gen aus dem Institut für angewandte Mathematik an der 
ETH. Zürich Nr. 3) 45 S. m. 8 Abb. und 3 Struktur- 
diagrammen. Basel 1952. Verlag Birkhäuser. Preis 
brosch. 5,70 Schw. Fr. 


Dr.-Ing. W.Matz, Aufgabensammlungzur 
Thermodynamik des Wärme- und 
Stoffaustauschesin der Verfahrens- 
teehnik. XII+ 138 S. m. 29 Abb. Darmstadt 1953. 
Verlag D. Steinkopf. Preis brosch. 16,— DM, geb. 
18,— DM. 


Dr.-Ing. N. Rosenauer, (Prof. a. d. Univ. von New 
South Wales) und A. H. Willis, B. Sc. (Eng.), Ph. D., 


Wh. So. (Prof. a.d. Univ. von New South Wales). 
Kinematicsof Mechanisms. XV + 395 S. mit 
377 Abb. Sidney 1953. Associated General Publi- 
cations. Preis gebunden 4/10/0 £ Aust. 


Dr. I. G. Petrovskij (o. Prof. a. d. Univ. Moskau). 
Vorlesungen über die Theorie der 
Integralgleiehungen. (Übersetzt von Dipl.- 
Math. R. Herschel.) 100S. Würzburg 1953, 
Physika-Verlag. Preis 7,80 DM. 


M.N. Butschin, Reibungsuntersuchun- 
gen an geschmierten Oberflächen 
bei tiefen Temperaturen (mit einem Vor- 
wort von Prof. Dr.-Ing. Dr.-Ing. E. h. E.Heide- 
broek Dresden). VIII + 64 S. mit 46 Abb. Berlin 
1953. Akademie-Verlag. Preis gebunden 7,50 DM. 
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NACHRICHTEN 


Stellenvermittlung zur Förderung des wissenschaftlichen 
Nachwuchses auf den Gebieten der Mathematik und 
Mechanik 


Wiederholte Anfragen stellensuchender Nachwuchs- 
kräfte der Mathematik und Mechanik, aber auch An- 
fragen aus Forschung, Industrie und Wirtschaft nach 
speziell qualifizierten Wissenschaftlern haben den Vor- 
standsrat der GaMM in seiner Aachener Sitzung vom 
21. April 1953 veranlaßt, geeignete Schritte zu ver- 
suchen, um dem hier zutage getretenen Bedürfnis abzu- 
helfen. Nicht zuletzt war dabei auch die Erwägung von 
Bedeutung, daß die Hochschulen den Studierenden in 
langjähriger Ausbildung ihre ganze Kraft widmen, 
ihrem weiteren DBerufsweg nach Abschluß der 
Hochschulbildung jedoch nicht immer die gebührende 
Aufmerksamkeit schenken. Insbesondere scheint es 
notwendig, sich der Situation der Diplom-Mathematiker 
zu widmen, deren Wirkungsbereich in der Praxis sich 
zwar ständig ausweitet, jedoch nach wie vor ent- 
scheidend von der Initiative einer jeden einzelnen in 
Industrie, Wirtschaft oder Verwaltung eintretenden 
Nachwuchskraft bestimmt wird und in seiner allge- 
meinen Anerkennung in der Öffentlichkeit noch man- 
ches zu wünschen übrig läßt. In noch erhöhtem Maße 
muß die Sorge der GaMM dem wissenschaftlichen Nach- 
wuchs auf den Hochschulen und in den Forschungs- 
instituten gelten. 

Aus diesen Erwägungen beraus hat der Vorstands- 
rat der GaMM die Geschäftsstelle der Gesellschaft (Frei- 
burg/Br. Institut für angewandte Mathematik der Uni- 
versität) beauftragt, eineStellenvermittlung für den jun- 
gen Nachwuchs einzurichten. Es wird dabei angenom- 
men, daß nach einer gewissen Anlaufzeit die erreichte 
Übersicht über Angebot und Nachfrage es gestatten 
wird, dem Nachwuchs den Weg zu einer der jeweiligen 
speziellen Eignung und den besonderen Wünschen ent- 
sprechenden Bestätigung zu ebnen und zu verkürzen, 
und daß andererseits eine zunebmende Aufklärung der 
Industrie und Wirtschaft über die Wirkungsmöglich- 
keiten der vorhandenen Nachwuchskräfte ihre Früchte 
tragen wird. Die neue Einrichtung soll aber auch und 
vor allem dazu beitragen, daß hochqualifizierte Kräfte, 
die an ihrer eigenen Hochschule keine freien Assistenten- 
stellen vorfinden, der Forschung durch Vermittlung an 
andere Hochschulen oder Forschungsinstitute erhalten 


bleiben. Sie könnte aber auch dazu beitragen, eine ge- 
wisse Beweglichkeit und den oft heilsamen „Klima- 
wechsel“ unter den Assistenten zu fördern. 

In Ausführung des genannten Beschlusses gibt die 
Geschäftsstelle der GaMM zunächst bekannt, daß 
stellensuchende Nachwuchskräfte mit abgeschlossener 
Hochschulbildung ihre Unterlagen (Lebenslauf mit 
Darstellung ibres wissenschaftlichen Werdegangs, ihrer 
besonderen wissenschaftlichen Betätigungen und Nei- 
gungen, Zeugnisabschriften sowie jede andere zweck- 
dienliche Angabe) bei der Geschäftsstelle der GaMM 
einreichen können. Die Geschäftsstelle der GaMM 
wird ihrerseits bestrebt sein, alle in Frage kommenden 
Institute und wissenschaftlichen Organisationen sowie 
geeignete Stellen der Industrie und Wirtschaft zu 
unterrichten. und sie auffordern, gegebenenfalls 
Stellenangebote mit Angabe der geforderten speziellen 
Qualifikationen an die GaMM gelangen zu lassen. 

Die Geschäftsstelle der GaMM sieht ihre Aufgabe 
darin, die erstrebten Verbindungen nach Maßgabe der 
vorbandenen Möglichkeiten herzustellen und Ver- 
handlungen anzubabnen. Sie verwendet dabei das ihr 
zur Verfügung gestellte Material, schaltet sich jedoch 
nicht in die Verhandlungen selbst ein und gibt nach 
keiner Seite eigene gutachtliche Äußerungen ab. Für 
die Bewerber besteht die Verpflichtung, ihre Bewer- 
bung bei der Geschäftsstelle in dem Augenblick schrift- 
lich zurückzuziehen, in dem sie eine Anstellung mit oder 
ohne Vermittlung der Geschäftsstelle gefunden haben, 
die ihren Wünschen entspricht. 

Die Bewerber um Stellungen werden gebeten, an die 
von der GaMM zu ihrer Förderung geschaffene Ein- 
richtung keine unangemessenen Erwartungen zu 
knüpfen und insbesondere zur Kenntnis zu nehmen, 
daß die „Geschäftsstelle der GaMM“ lediglich über 
wenige ehrenamtlich tätige Kräfte verfügt. Die Be- 
arbeitung von Fragen einer allgemeinen Berufsberatung 
kann unter diesen Umständen nicht übernommen, und 
der Schriftverkebr muß auf ein Minimum beschränkt 
werden. Von der Mitwirkung der durch dieses Unter- 
nehmen Angesprochenen wird es abhängen, ob der 
Versuch einer Stellenvermittlung zu dem gewünschten 
Erfolg führt oder nicht. 


Freiburg i. Br. H. Görtler 
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